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Summary. — 

Comprehensive description of the phenomenology of the 13(3 decay is given, with em- 
phasis on the nuclear physics aspects. After a brief review of the neutrino oscillation 
results and of motivation to test the lepton number conservation, the mechanism of 
the Qvf3f3 is discussed. Its relation to the lepton flavor violation involving charged 
leptons and its use as a diagnostic tool of the Qvl3(3 mechanism is described. Next 
the basic nuclear physics of both ,5/3-decay modes is presented, and the decay rate 
formulae derived. The nuclear physics methods used, the nuclear shell model and 
the quasiparticle random phase approximation, are described next, and the choice 
of input parameters is discussed in the following section. Finally, the numerical 
values of the nuclear matrix elements, and their uncertainty, are presented. In the 
appendix the relation of the search for the neutrino magnetic moment to the Dirac 
versus Majorana nature of neutrinos is described. 
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1. Introduction to /3/3 decay 

In the last decade neutrino oscillation experiments have convincingly and triumphantly 
shown that neutrinos have a finite mass and that the lepton flavor is not a conserved 
quantity. These results opened the door to what is often called the "Physics Beyond 
the Standard Model" . In other words, accommodating these findings into a consistent 
scenario requires generalization of the Standard Model of electroweak interactions that 
postulates that neutrinos are massless and that consequently lepton fiavor, and naturally, 
also the total lepton number, are conserved quantities. 

In oscillation experiments only the differences in squares of the neutrino masses, 
Am^ = |m| — TO^I, is measured, and the results do not depend on the charge conjugation 
properties of neutrinos, i.e. whether they are Dirac or Majorana fermions. Nevertheless, 
from oscillation experiments one can establish a lower limit on the absolute value of the 
neutrino mass scale, niscaie = -y/IAm^l. Thus, one or more neutrinos have a mass of at 
least ^|Am^j„| ~ 50 meV, and another one has mass of at least ^\ls.m1^i \ ~ 10 meV. 
In addition, an upper limit on the masses of all active neutrinos ~ 2 - 3 eV can be derived 
from the combination of analysis of the tritium beta-decay experiments and the neutrino 
oscillation experiments. Combining these constraints, masses of at least two (out of the 
total of three active) neutrinos are bracketted by 10 meV < < 2 - 3 eV. 

Thus, neutrino masses are six or more orders of magnitude smaller than the masses 
of the other fermions. Moreover, the pattern of masses, i.e. the mass ratios of neutrinos, 
is rather different (even though it remains largely unknown) than the pattern of masses 
of the up- or down- type quarks or charged leptons. All of these facts suggest that, 
perhaps, the origin of the neutrino mass is different than the origin (which is still not 
well understood) of the masses of the other fermions. 

The discoveries of neutrino oscillations, in turn, are causing a renaissance of enthusi- 
asm in the double beta decay community and a slew of new experiments that are expected 
to reach, within a near future, the sensitivity corresponding to the neutrino mass scale. 
Below I review the current status of the double beta decay and the effort devoted to 
reach the required sensitivity, as well as various issues in theory (or phenomenology) 
related to the relation of the Qv[3(3 decay rate to the absolute neutrino mass scale and to 
the general problem of the Lepton Number Violation (LNV). And, naturally, substantial 
emphasis is devoted to the nuclear structure issues. But before doing that I very briefiy 
summarize the achievements of the neutrino oscillation searches and the role that the 
search for the neutrinoless double beta decay plays in the elucidation of the pattern of 
neutrino masses and mixing. 

There is a consensus that the measurement of atmospheric neutrinos by the Su- 
perKamiokande collaboration [1^ can be only interpreted as a consequence of the nearly 
maximum mixing between z/^ and Vt neutrinos, with the corresponding mass squared 
difference lAm^t^l ~ 2.4 X lO-^eV^. This finding was confirmed by the K2K experiment 
[2] that uses accelerator beam pointing towards the SuperKamiokande detector 250 
km away as well as by the very recent result of the MINOS experiment located at the 
Sudan mine in Minnesota 735 km away from Fermilab|3^. Several large long-baseline 
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Fig. 1. - Schematic illustration of the decomposition of the neutrino mass eigenstates Vi in 
terms of the flavor eigenstates. The two hierarchies cannot be, at this time, distinguished. The 
small admixture of into vz is an upper limit, and the mass square of the neutrino the 



quantity mf , remains unknown. 



experiments are being build to further elucidate this discovery, and determine the corre- 
sponding parameters even more accurately. 

At the same time the " solar neutrino puzzle" , which has been with us for over thirty 
years since the pioneering chlorine experiment of Davis@], also reached the stage where 
the interpretation of the measurements in terms of oscillations between the t'e and some 
combination of active, i.e., v^^ and Vr neutrinos, is inescapable. In particular, the jux- 
taposition of the results of the SNO experiment 5 and SuperKamiokande |B] , together 
with the earlier solar neutrino flux determination in the gallium experiments [7l and, 
of course chlorine 4 , leads to that conclusion. The value of the corresponding oscillation 
parameters, however, remained uncertain, with several "solutions" possible, although the 
so-called Large Mixing Angle (LMA) solution with sin^ 20soi ~ 0.8 and Am^^j ^ lO^^'eV^ 
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was preferred. A decisive confirmation of the "solar" oscillations was provided by the 
nuclear reactor experiment KamLAND [S] [TOl ITT] that demonstrated that the flux of the 
reactor v,, is reduced and its spectrum distorted at the distance Lq 180 km from nu- 
clear reactors. The most recent KamLAND results [11 , combined with the existing solar 
neutrino data launched the era of precision neutrino measurements, with the correspond- 
ing parameters Am^^ = 7.59lQ;2i x lO^^eV^ and tan^ 6*12 = 0.47tQ;o5 determined with 
an unprecedented accuracy. Analysis of that experiment, moreover, clearly shows the 
oscillatory behavior of the detection probability as a function of Lq/Eu. That behavior 
can be traced in ref. jllj over two full periods. 

The pattern of neutrino mixing is further simplified by the constraint due to the 
Chooz and Palo Verde reactor neutrino experiments [121 [13] which lead to the conclusion 
that the third mixing angle, ^13, is small, sin^ 26*13 < 0.1. The two remaining possible 
neutrino mass patterns are illustrated in Figjl] 

As already stated, oscillation experiments cannot determine the absolute magnitude of 
the masses and, in particular, cannot at this stage separate two rather different scenarios, 
the hierarchical pattern of neutrino masses in which m ~ V Am? and the degenerate 
pattern in which m ^ V Am^. It is hoped that the search for the neutrinoless double 
beta decay, reviewed here, will help in foreseeable future in determining or at least 
narrowing down the absolute neutrino mass scale, and in deciding which of these two 
possibilities is applicable. 

Moreover, even more important is the fact that the oscillation results do not tell us 
anything about the properties of neutrinos under charge conjugation. While the charged 
leptons are Dirac particles, distinct from their antiparticles, neutrinos may be the ulti- 
mate neutral particles, as envisioned by Majorana, that are identical to their antiparticles. 
That fundamental distinction becomes important only for massive particles. Neutrinoless 
double beta decay proceeds only when neutrinos are massive Majorana particles, hence 
its observation would resolve the question. 

The argument for the "Majorana nature" of the neutrinos can be traced to the ob- 
servation by Weinberg [Tl] who pointed out almost thirty year ago that there exists only 
one lowest order (dimension 5, suppressed by only one inverse power of the corresponding 
high energy scale A) gauge-invariant operator given the content of the standard model 



where L is the lepton doublet, L'^ = L^C with C the charge conjugation operator, 
e = — ZT2, and H represent the Higgs boson. After the spontaneous symmetry breaking 
the Higgs acquires vacuum expectation value and the above operator represents the 
neutrino Majorana mass that violates the total lepton number conservation law by two 
units 



(1) 



£(5) =c(5)/A(L^ei/)(i/^eL) , 



(2) 



A 2 



(i/^i^) -f- h.c. , 



where v ^ 250 GeV and the neutrinos are naturally light because their mass is suppressed 
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by the large value of the new physics scale A in the denominator. 

The most popular explanation of the smallness of neutrino mass is the see-saw mecha- 
nism, which is also roughly thirty years old [15 . In it, the existence of heavy right-handed 
neutrinos Nn is postulated, and by diagonalizing the corresponding mass matrix one ar- 
rives at the formula 



where the Dirac mass mu is expected to be a typical charged fermion mass and Mm is 
the Majorana mass of the heavy neutrinos Nn. Again, the small mass of the standard 
neutrino is related to the large mass of the heavy right-handed partner. Requiring that 
TOy is of the order of 0.1 eV means that Mm (or A) is ^ 10^'*"^^ GeV, i.e. near the GUT 
scale. That makes this template scenario particularly attractive. 

Clearly, one cannot reach such high energy scale experimentally. But, these scenarios 
imply that neutrinos are Majorana particles, and consequently that the total lepton 
number should not be conserved. Hence the tests of the lepton number conservation 
acquires a fundamental importance. 

Double beta decay is a nuclear transition (Z, A) ^ {Z + 2, A) in which two 

neutrons bound in a nucleus are simultaneously transformed into two protons plus two 
electrons (and possibly other light neutral particles). This transition is possible and 
potentially observable because nuclei with even Z and N are more bound than the odd- 
odd nuclei with the same A ~ N + Z. Analogous transition of two protons into two 
neutrons are also, in principle, possible in several nuclei, but phase space considerations 
give preference to the former mode. 

An example is shown in Fig. [2] The situation shown there is not really exceptional. 
There are eleven analogous cases (candidate nuclei) with the Q- value (i.e. the kinetic 
energy available to leptons) in excess of 2 MeV. 

There are two basic modes of the f3/3 decay. In the two-neutrino mode (2vl3(3) there 
are 2 Vf. emitted together with the 2 e~. It is just an ordinary beta decay of two bound 
neutrons occurring simultaneously since the sequential decays are forbidden by the energy 
conservation law. For this mode, clearly, the lepton number is conserved and this mode of 
decay is allowed in the standard model of electroweak interaction. It has been repeatedly 
observed in a number of cases and proceeds with a typical half-life of ~ lO^^^^^years 
for the nuclei with Q-values above 2 MeV. In contrast, in the neutrinoless mode {Qv(3(3) 
only the 2e~ are emitted and nothing else. That mode clearly violates the law of lepton 
number conservation and is forbidded in the standard model. Hence, its observation 
would be a signal of a "new physics". 

The two modes of the /3/3 decay have some common and some distinct features. The 
common features are: 



• The leptons carry essentially all available energy. The nuclear recoil is negligible, 
Q/Arup < 1. 
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Fig. 2. - Atomic masses of the isotopes with A = 136. Nuclei "'^Xe, "®Ba and ""Ce are stable 
against the ordinary /3 decay; hence they exist in nature. However, energy conservation alone 
allows the transition ^^^Xe ^"^^Ba + 2e~ (+ possibly other neutral light particles) and the 
analogous decay of ^^®Ce with the positron emission. 



• The transition involves the 0+ ground state of the initial nucleus and (in almost all 
cases) the 0+ ground state of the final nucleus. In few cases the transition to an 
excited 0+ or 2+ state in the final nucleus is energetically possible, but suppressed 
by the smaller phase space available. (But the 2z//3/3 decay to the excited O"*" state 
has been observed in few cases.) 

• Both processes are of second order of weak interactions, ^ Gp, hence inherently 
slow. The phase space consideration alone (for the 2i^/3/? mode ~ Q^^ and for the 
OvPP mode ~ Q^) give preference to the 01^/3/3 which is, however, forbidden by the 
lepton number conservation. 

The distinct features are: 

• In the 2uPj3 mode the two neutrons undergoing the transition are uncorrelated (but 
decay simultaneously) while in the Qvj3l3 the two neutrons are correlated. 

• In the 2vl3(3 mode the sum electron kinetic energy Ti + T2 spectrum is continuous 
and peaked below Q/2. This is due to the electron masses and the Coulomb 
attraction. As Ti + T2 ^ Q the spectrum approaches zero approximately like 



• On the other hand in the Qv(i(3 mode the sum of the electron kinetic energies is 
fixed, T1 + T2 — Q, smeared only by the detector resolution. 
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Fig. 3. - Separating the Qh/fSf} mode from the 2uj3l3 by the shape of the sum electron spectrum, 
including the effect of the 2% resolution smearing. The figure is for the rate ratio 1/100 and 
the insert for 1/10*^. 



These last distinct features allow one to separate the two modes experimentally by 
measuring the sum energy of the emitted electrons with a good energy resolution, even 
if the decay rate for the Qv[3(3 mode is much smaller than for the 2vl3(3 mode. This 
is illustrated in Figjs] where the insert shows the situation for the rate ratio of 1 : 10^ 
corresponding to the most sensitive current experiments. 

Various aspects, both theoretical and experimental, of the /3/3 decay have been re- 
viewed many times. Here I quote just some of the review articles [TBI [Tfl [T51 [TOl EUl 1^ . 
earlier references can be found there. 

In this introductory section let me make only few general remarks. The existence of 
the decay would mean that on the elementary particle level a six fermion lepton 

number violating amplitude transforming two d quarks into two u quarks and two elec- 
trons is nonvanishing. As was first pointed out by Schechter and Valle [22] more than 
twenty years ago, this fact alone would guarantee that neutrinos are massive Majorana 
fermions (see Fig. |4]). This qualitative statement (or theorem), however, does not in 
general allow us to deduce the magnitude of the neutrino mass once the rate of the 0v(3p 
decay have been determined. It is important to stress, however, that quite generally an 
observation of any total lepton number violating process, not only of the 0i/(3/3 decay, 
would necessarily imply that neutrinos are massive Majorana fermions. 

There is no indication at the present time that neutrinos have nonstandard interac- 
tions, i.e. they seem to have only interactions carried by the W and Z bosons that are 
contained in the Standard Electroweak Model. All observed oscillation phenomena can 
be understood if one assumes that neutrinos interact exactly the way the Standard Model 
prescribes, but are massive fermions forcing a generalization of the model. If we accept 
this, but in addition assume that neutrinos are Majorana particles, we can in fact relate 
the QvPP decay rate to a quantity containing information about the absolute neutrino 
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Fig. 4. - By adding loops involving only standard weak interaction processes the /3/3 decay 
amplitude (the black box) implies the existence of the Majorana neutrino mass. 



mass. With these caveats that relation can be expressed as 

(4) -^^G°''(Q,Z)|M°''P(m^^)^ 

^1/2 

where G'^'^{Q,Z) is a phase space factor that depends on the transition Q value and 
through the Coulomb effect on the emitted electrons on the nuclear charge and that can 
be easily and accurately calculated (a complete list of the phase space factors G'^'^{Q, Z) 
and G^'^{Q,Z) can be found, e.g. in Ref. |23]), M°'' is the nuclear matrix element that 
can be evaluated in principle, although with a considerable uncertainty and is discussed 
in detail later, and finally the quantity {rUfsp) is the effective neutrino Majorana mass, 
representing the important particle physics ingredient of the process. 

In turn, the effective mass (w/j/s) is related to the mixing angles 6ij (or to the matrix 
elements |C/e,i| of the neutrino mixing matrix) that are determined or constrained by the 
oscillation experiments, to the absolute neutrino masses rui of the mass eigenstates Vi 
and to the totally unknown additional parameters, as fundamental as the mixing angles 
9ij, the so-called Majorana phases a{i), 

(5) {m(30) = |S,|[/e,pe'"Wm,| . 

Here Uei are the matrix elements of the first row of the neutrino mixing matrix. 

It is straightforward to use eq.(|5| and the known neutrino oscillation results in order 
to compare (mpp) with other neutrino mass related quantities. This is illustrated in 
FigjS] Traditionally such plot is made as in the left panel. However, the lightest neutrino 
mass rrimin is not an observable quantity. For that reason the other two panels show the 
relation of {mpp) to the sum of the neutrino masses M that is constrained and perhaps 
one day will be determined by the "observational cosmology" , and also to (to/j) that 
represent the parameter that can be determined or constrained in ordinary [3 decay. 



(6) 



Nuclear physics aspects of double beta decay 



9 



Several remarks are in order. First, the observation of the Qv(3[3 decay and determi- 
nation of (mpp), even when combined with the knowledge of M and/or (m^) does not 
allow, in general, to distinguish between the normal and inverted mass orderings. This 
is a consequence of the fact that the Majorana phases are unknown. In regions in Fig. [5] 
where the two hatched bands overlap it is clear that two solutions with the same {rriisis) 
and the same M (or the same {mp)) always exist and cannot be distinguished. 

On the other hand, obviously, if one can determine that (n^/3/3) > 0.1 eV we would 
conclude that the mass pattern is degenerate. And in the so far hypothetical case that one 
could show that {mpp) < 0.01 eV, but nonvanishing nevertheless, the normal hierarchy 
would be establishec (^) 

It is worthwhile noting that if the inverted mass ordering is realized in nature, (and 
neutrinos are Majorana particles) the quantity (w/j^s) is constrained from below by ~0.01 
eV. This is within the reach of the next generation of experiments. Also, at least in prin- 
ciple, in the case of the normal hierarchy while all neutrinos could be massive Majorana 
particles it would be still possible that {rri/sis) = 0. Such a situation, however, requires 
"fine tuning" or reflects a symmetry of some kind. For example, if 6*13 = the relation 
TO1/TO2 = tan^ 012 must be realized in order that (mpp) — 0. This implies, therefore, a 
definite relation between the neutrino masses and mixing angles. There is only one value 
of mi — 4.6 meV for which this condition is valid. 

Let us finally remark that the Ov/S/S decay is not the only LNV process for which 
important experimental constraints exist. Examples of the other LNV processes with 
important limits are 

fi^ + {Z, A) ^ + {Z - 2, A); exp. branching ratio < 10"^^ , 
/i^'^/i+TT^; exp. branching ratio < 3 x 10~^ , 
(7) Ve emission from the Sun; exp. branching ratio < 10^** . 

However, detailed analysis suggests that the study of the Qvf3(i decay is by far the most 
sensitive test of LNV. In simple terms, this is caused by the amount of tries one can 
make. A 100 kg Qvj3f3 decay source contains ~ 10^^ nuclei that can be observed for a 
long time (several years). This can be contrasted with the possibilities of first producing 
muons or kaons, and then searching for the unusual decay channels. The Fermilab 
accelerators, for example, produce ~ 10^*^ protons on target per year in their beams and 
thus correspondingly smaller numbers of muons or kaons. 

2. Mechanism of the Qvj3(3 decay 

It has been recognized long time ago that the relation between the Oi'/S/^-decay rate 
and the effective Majorana mass (mpp) is to some extent problematic. The rather con- 



(^) In that case also the {mp) in the right panel would not represent the quatity directly related 
to the ordinary (3 decay. There are no ideas, however, how to reach the corresponding sensitivity 
in ordinary (3 decay at the present time. 
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Fig. 5. - The left panel shows the dependence of (mpp) on the absolute mass of the lightest 
neutrino mass eigenstate rrimin, the middle one shows the relation between {mpp) and the sum 
of neutrino masses M = Em^ determined or costrained by the "observational cosmology", and 
the right one depicts the relation between (m/j^) and the effective mass (m/j) determined or 
contrained by the ordinary /3 decay. In all panels the width of the hatched area is due to 
the unknown Majorana phases and therefore irreducible. The solid lines indicate the allowed 
regions by taking into account the current uncertainties in the oscillation parameters; they will 
shrink as the accuracy improves. The two sets of curves correspond to the normal and inverted 
hierarchies, they merge above about {m/s/s) ~ 0.1 eV, where the degenerate mass pattern begins. 



servative assumption leading to eq.Q is that the only possible way the 0iyf3f3 decay can 
occur is through the exchange of a virtual light, but massive, Majorana neutrino be- 
tween the two nucleons undergoing the transition, and that these neutrinos interact by 
the standard left-handed weak currents. But that is not the only possible mechanism. 
LNV interactions involving so far unobserved much heavier TeV) particles can lead 
to a comparable OvPf] decay rate. Some of the possible mechanisms of the elementary 
dd ^ uu + e^e^ transition (the "black box in Fig. |4| are indicated in Fig. |6] Only 
the graph in the upper left corner would lead to eq. Q. Thus, in the absence of addi- 
tional information about the mechanism responsible for the Qu(3j3 decay one could not 
unambiguously infer the magnitude of {mpp) from the Qv P (S-decay rate. 

In general decay can be generated by (i) light massive Majorana neutrino ex- 

change or (ii) heavy particle exchange (see, e.g. Refs.[241 |2S]), resulting from LNV 
dynamics at some scale A above the electroweak one. The relative size of heavy {Ah) 
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Fig. 6. - All these symbolic Feynman graphs potentially contribute to the Oz^/3/3-decay amplitude 



versus light particle (Al) exchange contributions to the decay ampHtude can be crudely 
estimated as follows 

^ ^^T^T' ^ AT ^ K^mpp) ' 

where {mpp) is the effective neutrino Majorana mass, (fc^) ~ (fOO MeV)^ is the typ- 
ical light neutrino virtuahty, and A is the heavy scale relevant to the LNV dynamics. 
Therefore, Ar/Al - 0(1) for {mpp) - 0.1 - 0.5 eV and A ~ 1 TeV, and thus the LNV 
dynamics at the TeV scale leads to similar Oj^/3/3-decay rate as the exchange of light 
Majorana neutrinos with the effective mass {mpf}) ~ 0.1 — 0.5 eV. 

Obviously, the lifetime measurement by itself does not provide the means for de- 
termining the underlying mechanism. The spin-flip and non-flip exchange can be, in 
principle, distinguished by the measurement of the single-electron spectra or polariza- 
tion (see e.g. [17]). However, in most cases the mechanism of light Majorana neutrino 
exchange, and of heavy particle exchange cannot be separated by the observation of the 
emitted electrons. Thus one must look for other phenomenological consequences of the 
different mechanisms. Here I discuss the suggestion [28 that under natural assumptions 
the presence of low scale LNV interactions, and therefore the absence of proportionality 
between {mpp)'^ and the Oi//3/3-decay rate also affects muon lepton flavor violating (LFV) 
processes, and in particular enhances the ^ e conversion compared to the — > 67 
decay. 

The discussion is concerned mainly with the branching ratios B^^e-y = r(/^ ~* 
e7)/r^°'' and iJ^^e = Tconv/rcapt, where /i — ^ 67 is normalized to the standard muon 
decay rate r|°' = {Glml) / {192n^) , while fj, e conversion is normalized to the capture 
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rate Fcapt- The main diagnostic tool in our analysis is the ratio 



(9) 



and the relevance of our observation relies on the potential for LFV discovery in the 
forthcoming experiments MEG [53] 67) and MECO [50] (a* e conversion) (•^) 

At present, the most stringent limit on the branching ratio B^^^^ is [31 1.2 x 10^^^ 
and the MEG experiment aims at the sensitivity about two orders of magnitude better. 
For the muon conversion the best experimental limit |32l used gold nuclei and reached 
Bf_i^e < 8 X 10^^"^. The various proposals aim at reaching sensitivity of about 10^^^. 

It is useful to formulate the problem in terms of effective low energy interactions 
obtained after integrating out the heavy degrees of freedom that induce LNV and LFV 
dynamics. Thus, we will be dealing only with the Standard Model particles and all sym- 
metry relations will be obeyed. However, operators of dimension > 4 will be suppressed 
by l/A''"^, where A is the scale of new physics. If the scales for both LNV and LFV 
are well above the weak scale, then one would not expect to observe any signal in the 
forthcoming LFV experiments, nor would the effects of heavy particle exchange enter 
Oi//3/3 at an appreciable level. In this case, the only origin of a signal in at the 

level of prospective experimental sensitivity would be the exchange of a light Majorana 
neutrino, leading to eq.Q, and allowing one to extract {rapp) from the decay rate. 

In general, however, the two scales may be distinct, as in SUSY-GUT |33] or SUSY 
see-saw models. In these scenarios, both the Majorana neutrino mass as well as LFV 
effects are generated at the GUT scale. The effects of heavy Majorana neutrino exchange 
in QvPP are, thus, highly suppressed. In contrast, the effects of GUT-scale LFV are 
transmitted to the TeV-scale by a soft SUSY-breaking sector without mass suppression 
via renormalization group running of the high-scale LFV couplings. Consequently, such 
scenarios could lead to observable effects in the upcoming LFV experiments but with an 
0{a) suppression of the branching ratio -B^_>e relative to -B^^e7 due to the exchange of 
a virtual photon in the conversion process rather than the emission of a real one. 

As a specific example let us quote the SUSY SU(5) scenario where [37] 



{^) Even though the MECO experiment, that aimed at substantial incerase in sensitivity of 
the — > e conversion, was recently cancelled, proposals for experiments with similar sensitivity 
exist elsewhere. 



(10) 




(11) 




where the gaugino masses were neglected. 
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Another example is the evaluation of the ratio TZ — B^^e/ B^^e-y in the constrained 
seesaw minimal supersymmetric model [33] with very high scale LNV and TZ 1/200 for 
a variety of input parameters. There are, however, exceptions, like the recent evaluation 
of 7?. in a variety of SUSY SO(IO) models '35] with high scale LNV but with TZ as large 
as 0.3 in one case. 

The case where the scales of LNV and LFV are both relatively low TeV) is more 
subtle and requires more detailed analysis. This is the scenario which might lead to 
observable signals in LFV searches and at the same time generate ambiguities in inter- 
preting a positive signal in Qv(3j3. This is the case where one needs to develop some 
discriminating criteria. 

Denoting the new physics scale by A, one has a LNV effective lagrangian (dimension 
d = 9 operators) of the form 

(12) CouPiJ Oi^ qViq qT2qeTze'' , 

i 

where we have suppressed the flavor and Dirac structures (a complete list of the dimension 
nine operators Oi can be found in Ref. [25 J. 

For the LFV interactions (dimension d = 6 operators) , one has 

(13) /:lfv - E 

i 

and a complete operator basis can be found in Refs. [351 ISH] • The LFV operators relevant 
to our analysis are of the following type (along with their analogues with L <-> i?): 

Oil — f-iL i'jL ^kL ^-rnL 

(14) Otq = liTdj qVqq . 

Operators of the type Oa are typically generated at one-loop level, hence our choice 
to explicitly display the loop factor l/(47r)^. On the other hand, in a large class of 
models, operators of the type Oi or Oiq are generated by tree level exchange of heavy 
degrees of freedom. With the above choices, all non-zero Ci are nominally of the same 
size, typically the product of two Yukawa-like couplings or gauge couplings (times flavor 
mixing matrices). 

With the notation established above, the ratio TZ of the branching ratios /i — > e to 
/i — > e-t-7 can be written schematically as follows (neglecting flavor indices in the effective 
coupHngs and the term with L ^ R) [28j : 



TZ 



487r2 



Ai e^Co-L + (A2QL + Mciq) log 



* a2 

• 2 I 2 

" 1^0/' fir -4— A ^ / •/! - 1 UlU 

ml 
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Fig. 7. - Loop contributions to fi —> e conversion througli insertion of operators Oe or Otq, 
generating the large logarithm. 



In the above formula Ai^2,3.4 are numerical factors of 0(1), while the overall factor 
arises from phase space and overlap integrals of electron and muon wavefunctions in 
the nuclear field. For hght nuclei $ = {ZF^)/{g^ + 3g\) ~ 0(1) (^y, A are the vector and 
axial nucleon form factors at zero momentum transfer, while Fp is the nuclear form factor 
at ~ ^iT^fi EH])- The dots indicate subleading terms, not relevant for our discussion, 
such as loop-induced contributions to q and ciq that are analytic in external masses 
and momenta. In contrast the logarithmically-enhanced loop contribution given by the 
second term in the numerator of TZ plays an essential role. This term arises whenever 
the operators OiL,R and/or 0(q appear at tree-level in the effective theory and generate 
one-loop renormalization of Oiq _3S] (see Fig.[7]). 

The ingredients in eq. ([T5| lead to several observations: (i) In absence of tree-level 
ct,L and C£q, one obtains TZ ~ (<I> a)/(127r) ^ 10~^ — 10^^, due to gauge couphng and 
phase space suppression, (ii) When present, the logarithmically enhanced contributions, 
i.e. when either ql or c^q or both are nonvanishing, compensate for the gauge coupling 
and phase space suppression, leading to 7?. ~ 0(1). (iii) If present, the tree-level coupling 
c^q dominates the /i — > e rate leading to 7?. 3> 1 . 

Thus, we can formulate our main conclusions regarding the discriminating power of 
the ratio TZ: 

1. Observation of both the LFV muon processes — > e and ^ ^ with relative 
ratio TZ ~ 10~^ implies, under generic conditions, that To^pp ~ (mpp)'^. Hence the 
relation of the 0i^/3f3 lifetime to the absolute neutrino mass scale is straightforward. 

2. On the other hand, observation of LFV muon processes with relative ratio TZ 3> 
10^^ could signal non-trivial LNV dynamics at the TeV scale, whose effect on 0i^/3p 
has to be analyzed on a case by case basis. Therefore, in this scenario no definite 
conclusion can be drawn based on LFV rates. 

3. Non-observation of LFV in muon processes in forthcoming experiments would imply 
either that the scale of non-trivial LFV and LNV is above a few TeV, and thus 
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^Oi^f30 ~ ("1/3/3)^, or that any TeV-scale LNV is approximately flavor diagonal (this 
is an important caveat). 

The above statements are illustrated using two explicit cases [5S]: the minimal su- 
persymmetric standard model (MSSM) with R-parity violation (RPV-SUSY) and the 
Left-Right Symmetric Model (LRSM) . Limits on the rate of the 0^(3(3 decay were used 
in the past to constrain parameters of these two models [24] . 

RPV SUSY — If one does not impose R-parity conservation [R — (— 
the MSSM superpotential includes, in addition to the standard Yukawa terms, lepton 
and baryon number violating interactions, compactly written as (see, e.g. [40] ) 

Wrpv = KjkLiLjEl + \[^f.LiQjDl + X'-j^.UiDjD'j. 
(16) +i^[UHu , 

where L and Q represent lepton and quark doublet superfields, while E'^ , U'^, D'^ are 
lepton and quark singlet superfields. The simultaneous presence of A' and A" couplings 
would lead to an unacceptably large proton decay rate (for SUSY mass scale Asusy ~ 
TeV), so we focus on the case of A" — and set /z' = without loss of generality. In 
such case, lepton number is violated by the remaining terms in Wj^pv, leading to short 



distance contributions to 0vP(3 [see Fig. |8ja)], with typical coefficients [cf. eq. (12l] 

Ci nas X'lii Tia-i Af^ 



(17) A. 



7 ^ 7 



where as,Q!2 represent the strong and weak gauge coupling constants, respectively. The 
RPV interactions also lead to lepton number conserving but lepton flavor violating op- 



erators [see Fig.^b)], with coefficients [cf. eq. (13)] 
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(18) 

where the flavor combinations contributing to can be found in Ref. |41j . Hence, for 
generic flavor structure of the couplings A and A' the underlying LNV dynamics generate 
both short distance contributions to Qvfifi and LFV contributions that lead to 7?. ^ 10^^, 
Existing limits on rare processes strongly constrain combinations of RPV couplings, 
assuming Ksusy is between a few hundred GeV and ^ 1 TcV. Non-observation of LFV 
at future experiments MEG and MECO could be attributed either to a larger Ksusv 
(> few TeV) or to suppression of couplings that involve mixing among first and second 
generations. In the former scenario, the short distance contribution to Qv(3(3 does not 
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Fig. 8. - Gluino exchange contribution to (a), and typical tree-level contribution to Ot 

(&) in RPV SUSY. 



compete with the long distance one [see eq. ([s])], so that To,yf3p ~ {mpp)'^. On the other 
hand, there is an exception to this "diagnostic tool". If the A and A' matrices are nearly 
flavor diagonal, the exchange of superpartners may still make non-negligible contributions 
to QvjSP without enhancing the ratio TZ . 

LRSM — The LRSM provides a natural scenario for introducing non-sterile, right- 
handed neutrinos and Majorana masses [42 . The corresponding electroweak gauge group 
SU{2)lxSU{2)rxU{1)b-l, breaks down to 5?7(2)l x C/(1)y at the scale A > ©(TeV). 
The symmetry breaking is implemented through an extended Higgs sector, containing a 
bi-doublet $ and two triplets A^^j^, whose leptonic couplings generate both Majorana 
neutrino masses and LFV involving charged leptons: 



(19) C'^^'^~U^ [y%<i> + yp) L], 

- (Ll)^' Li - (Z^^ y% A« Lj, . 

Here $ = o'2^*f2, ^l,r = io-2^L.R, and leptons belong to two isospin doublets L]^ ^ = 
i'^L fl'^L i?,)- The gauge symmetry is broken through the VEVs (A^) — v^, (A") = 0, 
(<&) — diag(Ki,K2). After diagonalization of the lepton mass matrices, LFV arises from 
both non-diagonal gauge interactions and the Higgs Yukawa couplings. In particular, 
the Ai^ij-lepton interactions are not suppressed by lepton masses and have the structure 
C ^ Aj"'^^^/iy (1 ± j5)£j + h.c. The couplings hij are in general non-diagonal and 
related to the heavy neutrino mixing matrix j43j. 

Short distance contributions to Oi^/3/3 arise from the exchange of both heavy vs and 
Al,r (see Fig.[9])(a), with 

f^2Q-j Ci g2 1 ^ 92 ^ee 



A5 M^^ M,, ' M^^ Ml 



where g2 is the weak gauge coupling. LFV operators are also generated through non- 
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Fig. 9. - Typical doubly charged Higgs contribution to Oi^/3/3 (a) and to Oi (b) in the LRSM. 



diagonal gauge and Higgs vertices, with (see Fig.|9jb)) 



hi,,h*^ Ca_ {h''h)ef, ._ 



Note that the Yukawa interactions needed for the Majorana neutrino mass necessarily 
imply the presence of LNV and LFV couplings hij and the corresponding LFV operator 
coefficients ce, leading to 7?, ~ 0(1). Again, non-observation of LFV in the next gener- 
ation of experiments would typically push A into the multi-TeV range, thus implying a 
negligible short distance contribution to Qv(3(3. As with RPV-SUSY, this conclusion can 
be evaded by assuming a specific flavor structure, namely yM approximately diagonal or 
a nearly degenerate heavy neutrino spectrum. 

In both of these phenomenologically viable models that incorporate LNV and LFV at 
low scale (~ TeV), one finds TZ ^ 10~^ p8l|4TJ|43]. It is likely that the basic mechanism 
at work in these illustrative cases is generic: low scale LNV interactions (AL = ±1 and/or 
AL = ±2), which in general contribute to QvfijS, also generate sizable contributions to 
/i — *■ e conversion, thus enhancing this process over /i — > 67. 

In conclusion of this section, the above considerations suggest that the ratio TZ = 
Bfi^e/ Bfj^^e-y of muou LFV processes will provide important insight about the mech- 
anism of neutrinoless double beta decay and the use of this process to determine the 
absolute scale of neutrino mass. Assuming observation of LFV processes in forthcoming 
experiments, if 7?. ~ 10~^ the mechanism of 0^/3/3 is light Majorana neutrino exchange 
and, therefore, I/T1/2 ^ {mpf})"^; on the other hand, liTZ^ 10^^, there might be TeV 
scale LNV dynamics, and no definite conclusion on the mechanism of Qu(3j3 decay can be 
drawn based only on LFV processes. 



3. Overview of the experimental status of the search for /3/3 decay 

Before embarking on the discussion of the nuclear structure aspects of the (3(3 decay 
let us briefly describe the experimental status of the field (more detailed information on 
this topics are in the lectures by A. Giuliani). The topic has a venerable history. The 
rate of the 2v(3f3 decay was first estimated by Maria Goeppert-Meyer already in 1937 
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in her thesis work suggested by E. Wigner, basically correctly. Yet, first experimental 
observation in a laboratory experiment was achieved only in 1987, fifty years later |44j . 
(Note that this is not really exceptional in neutrino physics. It took more than twenty 
years since the original suggestion of Pauli to show that neutrinos are real particles in 
the pioneering experiment by Raines and Cowan. And it took another almost fifty years 
since that time to show that neutrinos are massive fermions.) Why it took so long in 
the case of the (3P decay? As pointed out above, the typical half-life of the 2i//3/? decay 
is ~ 10^° years. Yet, its "signature" is very similar to natural radioactivity, present to 
some extent everywhere, and governed by the half-life of ~ 10^" years or much less for 
most of the man-made or cosmogenic radioactivities. So, background suppression is the 
main problem to overcome when one wants to study either of the /3/3 decay modes. 

During the last two decades the 2iyf3f3 decay has been observed in "live" laboratory 
experiments in many nuclei, often by different groups and using different methods. That 
shows not only the ingenuity of the experimentalists who were able to overcome the 
background nemesis, but makes it possible at the same time to extract the corresponding 
21^ nuclear matrix element from the measured decay rate. In the 2i/ mode the half-life is 
given by 



(22) l/ri/2 = G2''(Q,Z)|M2''|2 , 

where G^'^{Q, Z) is an easily and accurately calculable phase space factor. 

The resulting nuclear matrix elements M^"^, which have the dimension energy"^, are 
plotted in Fig 10 Note the pronounced shell dependence; the matrix element for ^°°Mo 



is almost ten times larger than the one for ^■^"Te. Evaluation of these matrix elements, 
to be discussed below, is an important test for the nuclear theory models that aim at the 
determination of the analogous but different quantities for the Ov ncutrinoless mode. 

The challenge of detecting the Oi^/3/3 decay is, at first blush, easier. Unlike the con- 
tinuous 21/(3(3 decay spectrum with a broad maximum at rather low energy where the 
background suppression is harder, the Qv(3(i decay spectrum is sharply peaked at the 
known Q value (see Figjs]), at energies that are not immune to the background, but a bit 
less difficult to manage. However, as also indicated in FigjS] to obtain interesting results 
at the present time means to reach sensitivity to the Qv half-lives that are ~ 10^ times 
longer than the 2v decay half-life of the same nucleus. 



The historical lessons are illustrated in Fig 11 where the past limits on the Oi^/3/3 decay 
half- lives of various candidate nuclei are translated using eq.Q into the limits on the 
effective mass {mpp). When plotted in the semi-log plot this figure represent "Moore's 
law" of double beta decay, and indicates that, provided that the past trend will continue, 
the mass scale corresponding to Arn^j„ will be reached in about 7 years. This is also 
the time scale of significant experiments these days. Note that the figure was made using 
some assumed values of the corresponding nuclear matrix elements, without including 
their uncertainty. For such illustrative purposes they are, naturally, irrelevant. 

The past search for the neutrinoless double beta decay, illustrated in Fig{TT] was 
driven by the then current technology and the resources of the individual experiments. 
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Fig. 10. - Nuclear matrix elements for the 2vPI3 decay extracted from the measured half-lives. 

The goal has been simply to reach sensitivity to longer and longer half-lives. The situation 
is different, however, now. The experimentalists at the present time can, and do, use 
the knowledge summarized in Fig. [5] to gauge the aim of their proposals. Based on 
that figure, the range of the mass parameter {mpp) can be divided into three regions of 
interest. 

• The degenerate mass region where all 3> y^Am^^. In that region {mpp) > 
0.1 eV, corresponding crudely to the Oiy half- lives of 10^6-27 years. To explore it 
(in a realistic time frame), ^ 100 kg of the decaying nucleus is needed. Several 
experiments aiming at such sensitivity are being built and should run very soon 
and give results within the next ~ 3 years. Moreover, this mass region (or a 
substantial part of it) will be explored, in a similar time frame, by the study of 
ordinary (3 decay (in particular of tritium, see the lectures by C. Weinheimer) and 
by the observational cosmology (see the lectures by S. Pastor). These techniques 
are independent on the Majorana nature of neutrinos. It is easy, but perhaps 
premature, to envision various possible scenarios depending on the possible outcome 
of these measurements. 

• The so-called inverted hierarchy mass region where 20 < (mpp) < 100 meV and the 
Oi^PP half-lives are about 10^'^"^^ years. (The name is to some extent a misnomer. 
In that interval one could encounter not only the inverted hierarchy but also a quasi- 
degenerate but normal neutrino mass ordering. SuccessfuU observation of the Qu[3j3 
decay will not be able to distinguish these possibilities, as I argued above. This is 
so not only due to the anticipated experimental accuracy, but more fundamentally 
due to the unknown Majorana phases.) To explore this mass region, ^ ton size 
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sources would be required. Proposals for the corresponding experiments exist, but 
none has been funded as yet, and presumably the real work will begin depending 
on the experience with the various ~ 100 kg size sources. Timeline for exploring 
this mass region is ~ 10 years. 

• Normal mass hierarchy region where {mpp) < 10-20 meV. To explore this mass 
region, ^ 100 ton sources would be required. There are no realistic proposals for 
experiments of this size at present. 

Over the last two decades, the methodology for double beta decay experiments has 
improved considerably. Larger amounts of high-purity enriched parent isotopes, com- 
bined with careful selection of all surrounding materials and using deep-underground 
sites have lowered backgrounds and increased sensitivity. The most sensitive experi- 
ments to date use ^^Ge, ^o^Mo, "^Cd, "Oxe, and ^^^Xe. For ^*^Ge the hfetime limit 
reached impressive values exceeding lO^^yearsjlSl US]. The experimental lifetime limits 
have been interpreted to yield effective neutrino mass limits typically a few eV and in 
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^^Ge as low as 0.3 - 1.0 eV (the spread reflects an estimate of the uncertainty in the 
nuclear matrix elements). Similar sensitivity to the neutrino mass was reached recently 
also in the CUORICINO experiment with "Oxe g?]. 

While all these experiments reported lower limits on the 0^/3/3-decay half-lives, a sub- 
set of members of the Heidelberg-Moscow collaboration |48] reanalyzed the data (and 
used additional information, e.g. the pulse-shape analysis and a different algorithm in 
the peak search) and claimed to observe a positive signal corresponding, in the latest 
publication, to the half-life T1/2 = 2.23^q32 ^ 10^^ years. That report has been followed 
by a lively discussion. Clearly, such an extraordinary claim with its profound implica- 
tions, requires extraordinary evidence. It is fair to say that a confirmation, both for the 
same ^^Ge parent nucleus, and better yet also in another nucleus with a different Q value, 
would be required for a consensus. In any case, if that claim is eventually confirmed (and 
the mechanism of the Oi^/3/3 decay determined, i.e. the validity of eq. Q assured), the 
degenerate mass scenario will be implicated, and eventual positive signal in the analysis 
of the tritium f3 decay and/or observational cosmology should be forthcoming. For the 
neutrinoless Pf3 decay the next generation of experiments, which will use ~ 100 kg of 
decaying isotopes will, among other things, test this recent claim. 



4. Basic nuclear physics of /3/3 decay 

Whether a nucleus is stable or undergoes weak decay is determined by the dependence 
of the atomic mass Ma of the isotope {Z, A) on the nuclear charge Z. This functional 
dependence near its minimum can be approximated by a parabola 



(A/2 — Z)^ 

(23) Ma{Z,A) = const -f 26,^„ ^42 + ^Couij^ + m^Z 



where the symmetry energy coefficient is hsym 50 MeV and the Coulomb energy 
coefficient is hcoui ~ 0.7 MeV. The m^Z term represents the mass of the bound electrons; 
their binding energy, for our purposes, is small enough to be neglected. The last term 
(5, decisive for the application to the /3/3 decay, describes nuclear pairing, the increase 
in binding as pairs of like nucleons couple to angular momentum zero. It is a small 
correction term and is given in a crude approximation by (5 ^ ± Ylj A^l"^ MeV for odd 
N and odd Z, or even N and even Z, respectively, while (5 for odd A. Thus, for odd 
A nuclei, typically only one isotope is stable; nuclei with charge Z smaller than the stable 
nucleus decay by electron emission, while those with larger Z decay by electron capture 
or positron emission or by both these modes simultaneously. For even A nuclei the 
situation is different. Due to the pairing term 5, the even-even nuclei form one parabola 
while the odd-odd nuclei form another one, at larger mass, as shown in Fig. [2] using A 
= 136 as an example. Consequently, in a typical case there exist two (or three as in Fig. 
[2]) even-even nuclei for a given A which are stable against both electron and positron (or 
EC) decays. As these nuclei usually do not have the same mass, the heavier may decay 
into the lighter through a second order weak process in which the nuclear charge changes 



by two units. This process is double beta decay. In Fig. 12 all nuclei that are practical 
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Q(MeV) Abuncl.(7o) 



48Ca^48Tj 


4.271 


0.187 


76ge -»^«'Se 


2.040 


7.8 




2.995 


9.2 




3.350 


2.8 




3.034 


9.6 


iiopd^iiocd 


2.013 


11.8 




2.802 


7.5 


124Sn-.124Te 


2.228 


5.64 


i30Te^i30Xe 


2.533 


34.5 




2.479 


8.9 




3.367 


5.6 



Fig. 12. - Candidate nuclei for (3(3 decay with 2e emission and with the Q value >2 MeV. The 
corresponding abundances are also shown. 

candidates for the search for the decay are listed. All of them exist in nature since 
their lifetime is longer than the age of the solar system. However, with a single exception 
(^•^"Te), all of them are relatively rare so that a large scale experiments also require a 
large scale and costly (and sometimes technically difficult) isotope enrichment. 

Double beta decay, therefore, proceeds between two even-even nuclei. All ground 
states of even-even nuclei have spin and parity 0+ and thus transitions 0+ — > 0+ are 
expected in all cases. Occasionally, population of the low-lying excited states of the 
daughter nucleus is energetically possible, giving rise to 0+ 2+ transitions or to tran- 
sitions to the excited 0+ states. 

Double beta decay with the electron emission (both the 2i/ and Oi^ modes) in which 
the nuclear charge increases by two units is subject to the obvious condition 

(24) Ma{Z, A) >Ma{Z + 2, A) , while Ma{Z, A) < Ma{Z +l,A) , 

where Ma is the atomic mass, with the supplementary practical requirement that single 
beta decay is absent, or that it is so much hindered (e.g., by the angular momentum 
selection rules) that it does not compete with double beta decay. 

Double beta decay with the positron emission and/or electron capture (both the 2i/ 
and Oz/ modes) in which the nuclear charge decreases by two units can proceed in three 
different ways: 

• Two positron emission when Ma{Z, A) > Ma{Z — 2, ^4) -h 4me , 

• One positron emission and one electron capture when Ma{Z, A) > Ma{Z — 2,A) + 
2me + Be , and 
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• Two electron captures (only the two neutrino mode, see below for the comment on 
the Oiy mode) when Ma{Z, A) > Ma{Z -2,A) + 8^(1) + B^{2) , 

where i?e is the positive binding energy of the captured electron. As mentioned above, 
a complete list of all candidate nuclei (except the double electron capture) with the 
corresponding phase space factors for both /3/3-decay modes is given in Ref . |53] . 

Since the relevant masses are atomic masses, the processes with emission of positrons 
have reduced phase space (terms with rUe above). While from the point of view of ex- 
perimental observation the positron emission seems advantageous (possibility to observe 
the annihilation radiation), the reduction of phase space means that the corresponding 
lifetimes are quite long. In fact, not a single such decay has been observed so far. 

The two electron capture decay without neutrino emission requires a special comment. 
Clearly, when the initial and final states have different energies, the process cannot pro- 
ceed since energy is not conserved. The radiative process, with bremsstrahlung photon 
emission, however, can proceed and its rate, unlike all the other neutrinoless processes, 
increases with decreasing Q value [49]. (However, the estimated decay rates are quite 
small and lifetimes long) . In the extreme case of essentially perfect degeneracy, a reso- 
nance enhancement can occur |50j . 

The case of resonance, though probably unrealistic, perhaps deserves some explana- 
tions. The initial state is the atom {Z,A), stable against ordinary (3 decay. The final 
state is the ion {Z — 2, A) with electron vacancies H, H' and, in general, with the nucleus 
in some excited state of energy E* . The resonance occurs if the final energy 

(25) E^E* +Eh + Eh' 

is close to the decay Q value, i.e. the difference of the initial and final atomic masses, and 
a perfect resonance occurs when Q — E is less than the width of the final state which is 
dominated by the electron hole widths Th,^h'- The decay rate neax resonance is given 
by the Breit-Wigner type formula 

(26) 1 - r 

where AM is the matrix element of weak interaction between the two degenerate atomic 
states. 

The states of definite energy, the eigenstates of the total hamiltonian, are superposi- 
tions of the initial and final states, mixed by AM. But in reality, the initial state is pure, 
and not a state of definite energy, since the final state decays essentially immediately. 

The mixing matrix element is |50j 

(27) AM ^ ^^7j'^^ (m,,)|^(0)pglM°- , 

where ip{0) is the amplitude at the origin of the wave function of the captured electrons 
and M"'' is the nuclear matrix element discussed later. Clearly, if the resonance can be 
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approached, the decay rate would be enhanced by the factor 4/r compared toT /{E—Q)^, 
where the width F is typically tens of eV. Estimates suggest that in such a case the decay 
lifetime for {rngrj) ^ 1 cV could be of the order of 10^"*"^^ years, competitive to the rate 
with 2e^ emission. However, chances of finding a case of a perfect (eV size) resonance 
when E is of order of MeV are very unlikely. 



5. — Decay rate formulae 

5'1. 2v decay. - Even though the 2i//3^-decay mode is unrelated to the fundamental 
particle physics issues, it is worthwhile to discuss it in some detail. This is so because it is 

the mode that is actually seen experimentally; it is also the inevitable background for the 
Oi'/J/J-decay mode. In nuclear structure theory the corresponding rate is used as a test 
of the adequacy of the corresponding nuclear models. And various auxiliary experiments 
can be performed to facilitate the evaluation of the M'^'^ nuclear matrix elements. 

The derivation of the 2i//3/3-decay rate formula is analogous to the treatment of ordi- 
nary beta decay. It begins with the Fermi golden rule for second order weak decay 



(28) ^ = 2Tr6{Eo-EfEf) 

where the sum over m includes all relevant virtual states in the intermediate odd-odd 

nucleus and f3 labels the different Dirac structures of the weak interaction Hamiltonian. 
Next we talce into account that the weak Hamiltonian is the product of the nuclear and 
lepton currents; the corresponding formula will include the summation over the indeces 
of the emitted leptons. Then the summation over the lepton polarizations is performed, 
taking into account the indistinguishability of the final lepton pairs. Because we are 
interested in the rate formula, we neglect terms linear in and that disappear after 
integration over angles. (The angular distribution of the electrons is of the form (1 - 
/3i • 02) for the 0+ 0+ transitions, where /Si is the velocity of the electron i). 
The energy denominators are of the following form 



-'m,/3 



{f\Hp\m}{m\HP\i) 



Ei — E„ 



E, 



(29) Km{Mm) = -= I ± ^ I , 

— Pvl + J^el -C/m — J^i + Pv2 + ■C'e2 

where the -|- sign belongs to and the — sign to and 

^ ^ m( m) _ ^ . _|_ -|- Ee\ — Ei + Pi,\ + Ee2 

The energy denominators in the factors K, M, L, N contain contributions of the nuclear 
energies E^ — Ei, as well as the lepton energies Ee+Pv ■ When calculating the 0+ — > 0+ 

transitions, it is generally a very good approximation to replace these lepton energies by 
the corresponding average value, i.e., Ee + p^ ^ -E'o/2, where the Eq = Mi — Mf is the 
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total decay energy including electron masses. In that case M„i = Nn = and 

^^^^ " ^™ ^ Era-E, + Eo/2 " ~ (M, + Mf)/2 ' 

The lepton momenta, for both electrons and neutrinos are all q < Q and thus qR ^ 1, 
where R is the nuclear radius. Hence the so-called long wavelength approximation is valid 



and the rate formula (with eq. (31 1 ) separates into a product of the nuclear and lepton 



parts, where the lepton part contains just the phase space integral 

(32) / F{Z,Eel)pelEeldE,i F{Z,E,2)Pe2Ee2dE,2{E^- E,i- Ee2f /^^ . 

where the integration over the neutrino momentum was already performed. The single 
electron spectrum is obtained by performing integration over dEf,2, while the spectrum 
of summed electron energies is obtained by changing the variables to E^i + Ee2 and 
Eel — Ee2 and performing the integration over the second variable. If an accurate result 
is required, the relativistic form of the function F(Z, E) must be used and numerical 
evaluation is necessary. 

For a qualitative and intuitive picture, one can use the simplified nonrelativistic 
Coulomb expression, so-called Primakoff-Rosen approximation [S2] 

(33) F(Z,£;)=^ 2^^" 



p 1 — 



This approximation allows us to perform the required integrals analytically. 

For example, the sum electron spectrum, which is of primary interest from the exper- 
imental point of view is then independent of Z, 

(34) S~^'-(^»-^'->=(l + ^*-+T + f + W 

where K is the sum of the kinetic energies of both electrons, in units of electron mass. 
The Coulomb effects result in shifting the maximum of dN /dK towards lower energy and 
in making the approach of dN/dK to zero when K ^ Tq steeper. 

The nuclear structure information is contained in the nuclear matrix element; only 
the Gamow- Teller ar part contributes in the long wavelength approximation 



(35) M 



2._^ (0;i'?.T+|™)(m|<7fcT+|0^ 



E„, - {A'U + Mf)/2 



The individual terms in the eq. ( 35 1 have a well defined meaning, in particular for 
the most relevant ground state to ground state transitions. The terms {'m\akT^\Qll) 
represent the [3~ strength in the initial nucleus and can be explored in the nucleon 
exchange reactions such as (p, n) and i^He, t). On the other hand the terms {(y^\(JiT^\'m) 
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represent the /S"*" strength in the final nucleus and can be explored in the nucleon charge 
exchange reactions such as {n,p) and (d,^ He). In this way one can (up to the sign) 
explore the contribution of several low lying states to the M'^'^ matrix element. 

It turns out that in several nuclei the lowest (or few lowest) 1+ states give a dominant 
contribution to M^"^. (This is so-called "single state dominance".) In those cases the 
above mentioned experiments allow one to determine the M'^^ indirectly, independently 
of the actual 2v[3l3 decay. Such data are, naturally, a valuable testing ground of nuclear 
theory. On the other hand, it is not a priori clear and easy to decide in which nuclei the 



sum over the l"*" states in eq. (35) converges very fast and in which nuclei many states 
contribute. We will return to this issue later. 

Since, as stated in Section 3, the half-lives of many 2v(}(3 decays was experimentally 
determined. One can then extract the values of the nuclear matrix elements M'^'^ using 



eq. (35 1 . They are depicted in Fig. 10 For completeness, we add a table of the most 



recent half-life measurements of the 2z//3/3 decay. 

Table I. - Summary of experimentally measured 2v[3j3 half-lives and matrix elements, mostly 
from the NEMO experiment J5ij / (^^^ Xe is an important exception where a limit is quoted). 



Isotope 


T?;2 (y) 


a4't (Mcv^^) 


*«Ca 


(3.9 ±0.7 ±0.6) X 10^^ 


0.05 ±0.01 




(1.7 ±0.2) X 10^^ 


0.13 ±0.01 


»2Se 


(9.6 ± 0.3 ± 1.0) X 10^^ 


0.10 ±0.01 


36Zr 


(2.0 ±0.3 ±0.2) X 10^^ 


0.12 ±0.02 


i"°Mo 


(7.11 ±0.02 ±0.54) X 10^** 


0.23 ±0.01 




(2.8 ±0.1 ±0.3) X lO^'' 


0.13 ±0.01 


128rpg{l) 


(2.0 ±0.1) X 10^* 


0.05 ± 0.005 




(7.6 ± 1.5 ±0.8) X 10^'' 


0.032 ±0.003 




> 1.0 X 10^2 (90% CL) 


<0.01 




(9.2 ± 0.25 ± 0.73) X 10^® 


0.06 ±0.003 


238^(2) 


(2.0 ±0.6) X 10^^ 


0.05 ±0.01 



'^'deduced from the geochemically determined half-life ratio ^^^Te/^^^Te 
'^'radiochemical result for all decay modes 



5'2. Qv rate. - We shall now indicate the derivation of the electron spectra and decay 
rates associated with the nonvanishing value of . The decay rate is of the general form 

(36) Wo^ = 27rS,p„|i?o^|2,5(i;el ± Ee2 +Ef- M,)d^Peld^Pe2 , 

where Ef is the energy of the final nucleus and Rq^, is the transition amplitude including 
both the lepton and nuclear parts. 
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The lepton part of the ampHtude is written as a product of two left-handed currents 

(37) e(a;)7p^(l - J5)'^j{x)e{y)-/„^{1 ~ j5)My) , 

where, Vj, represent neutrino mass eigenstates j and k, and there is a contraction over 
the two neutrino operators. The contraction above is allowed only if the neutrinos are 
Majorana particles. 

After substitution for the neutrino propagator and integration over the virtual neu- 
trino momentum, the lepton amplitude acquires the form 

(38) -iS,k J g2,^2 e(^)7p^(l-75)(g^7M + "^j-)2(^-^5)7>Te^(?/) • 
From the commutation properties of the gamma matrices it follows that 

(39) ^(l-75)('Z^7M + ™.)^(l-75)=m,i(l-75) • 

Thus the decay amplitude for purely left-handed lepton currents is proportional to the 
neutrino Majorana mass rrij. 

Integration over the virtual neutrino energy leads to the replacement of the propagator 
{q^ — m|)~^ by the residue n jujj with lOj = (q ^ + m|)^/^. For the remaining integration 
over the space part dq we have to consider, besides this denominator Wj, the energy 
denominators of the second order perturbation expression. Denoting 



(40) = S„ -E, + E, 



e I 



we find that integration over dq leads to an expression representing the effect of the 
neutrino propagation between the two nucleons. This expression has the form of a "neu- 
trino potential" and appears in the corresponding nuclear matrix elements, introducing 
dependence of the transition operator on the coordinates of the two nucleons, as well as 
a weak dependence on the excitation energy — Ei of the virtual state in the odd-odd 
intermediate nucleus. 

There are several neutrino potentials as explained in the next subsection. The main 
one is of the form 

(41) n(r.^ E,.) = « / = ^ r . 

27r^5^ J uj UJ + A,n TTrg% Jq u;{u; + Am) 

Here we added the nuclear radius R — 1.2A^/'^ fm as an auxiliary factor so that H 
becomes dimensionless. A corresponding 1/R^ compensates for this auxiliary quantity 
in the phase space formula. Note that a consistency between the definitions of the nuclear 



radius R is required. ) The first factor to in the denominator of eq. (411 is the residue 



while the factor uj + Am is the energy denominator of perturbation theory. To obtain the 
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final result one has to treat properly the antisymmetry between the identical outgoing 
electrons (see [55]). 

The momentum of the virtual neutrino is determined by the uncertainty relation 
q ^ l/r, where r < i? is a typical spacing between two nucleons. We will show later that 
in fact the relevant values of r are only r < 2-3 fm, so that the momentum transfer q ~ 
100-200 MeV. For the light neutrinos the neutrino mass rrij can then be safely neglected 
in the potential H{r). (Obviously, for heavy neutrinos, with masses Mj ^ 1 GeV a 
different procedure is necessary.) Also, given the large value of q the dependence on the 
difference of nuclear energies Em — Ei is expected to be rather weak and the summation 
of the intermediate states can be performed in closure for convenience. (We will discuss 
the validity of that approximation later.) 

Altogether, we can rewrite the expression for the neutrino potential as 

(42) H{r) = -<^{ur) , 

r 

where <I>(u;r) < 1 is a relatively slowly varying function of r. From that it follows that 
a typical value of H{r) is larger than unity, but less than 5-10. In the next subsection 
we present the exact expressions for the neutrino potentials and the Qvfip transition 
operator for the most interesting case of small but finite neutrino Majorana masses nij. 
For now we use the relation (|4| and evaluate the phase space function G""^ 

(43) G°'^(Q, Z)^ j F{Z, E,^)F{Z, E,2)PeiPe2EeiE,25{Eo - -Bei - E,2)dE,idE,2 ■ 
The constant factor in front of this expression is 



(44) {GFCosOcgA) 



ficV I 1 



R J ?iln(2)327r5 ' 



and the values of 1 /G'^^ listed in Ref. [2F are in years provided that the neutrino masses 
are in eV. 



Again, in the Primakoff-Rosen approximation, eq. (33 1, G^"^ is independent of Z and 
(only the Eq dependence is shown) 

(45) a',„^{§-'-f.E„-l 

where Eq is expressed in units of electron mass. Each of the two electrons observed 
separately will have energy spectrum determined by the phase space integral. In the 
Primakoff-Rosen approximation its shape is 

(46) ^^(Te+l)2(^o-Te+l)^ 
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where again the kinetic energies are in the units of electron mass. 

It is of interest to contrast the Oi/ and 2i/ decay modes from the point of view of 
the phase space integrals. The Oi^ mode has the advantage of the two-lepton final state, 
with the characteristic Eq dependence compared to the four-lepton final state with E^^ 
dependence for the 21^ mode. In addition, the large average momentum of the virtual 
neutrino, compared with the typical nuclear excitation energy also makes the Ov decay 
faster. Thus, if {mpp) were to be of the order of me, the Qv decay would be ~ 10^ times 
faster than the 2v decay. It is this phase space advantage which makes the Qvf3(3 decay 
a sensitive probe for Majorana neutrino mass. 

Finally, let us remark that it is possible that, in addition to the two electrons, a light 
boson, the so called majoron, is emitted. This transition would have three-body phase 
space, giving rise to a continuous spectrum peaked at approximately three quarters of 
the decay energy Tq. We shall not discuss this topic here, but refer to Doi et al. [521 for 
a discussion of this issue. 

Also, since we concentrate on the mechanism involving the exchange of light Majorana 
neutrinos, we will not discuss in any detail the case of hypothetical right-handed currents. 

5'3. Exact expressions for the transition operator. - After the qualitative discussion in 
the preceding subsection, here we derive the exact expressions for the transition operator. 
The hadronic current, expressed in terms of nucleon fields ^, is 



where mp is the nucleon mass and is the momentum transfer, i.e. the momentum of 
the virtual neutrino. Since in the QvjSP decay ^ Qq we take ~ ~q^- 

For the vector and axial vector form factors we adopt the usual dipole approximation 



with gv = 1, gA = 1-254, My = 850 MeV, and Ma = 1086 MeV. These form factors 
are a consequence of the composite nature of nucleons. With high momentum transfer 
the "elastic" transitions, in which a nucleon remains nucleon and no other hadrons are 
produced, is reduced. 

We use the usual form for the weak magnetism, and the Goldberger-Treiman relation 
for the induced pseudoscalar term: 



(47) JP^ = ^T+ gy{q^)Y + igM{q^) q^-gA{q^)Yj5-9p{q^)q'l5 



(48) 



9v{q^) = 9v/{l + qVM^)^, gA{q^) ^ Qa/ + q^ I Mlf , 



(49) 



gniq'^) = {f^p- fJ.n)9v{q^) , gp{q^) = 2mpgA{q^)/{q^ + ml) . 



Reducing the nucleon current to the non-relativistic form yields (see Ref . |54] ) : 



A 



(50) 




n=l 



k 
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where J°{q^) = gv{q^) and 



(51) 



Jn{q ) ^ gAiiq \- 9A{q )o - gp{q )- 



2m„ 



2m„ 



r„ is the coordinate of the nth nucleon, k = 1,2,3, and gP'"' is the metric tensor. 

This allows us to derive the effective two-body transition operator in the momentum 
representation: 



(52) n 



- {—hp + hGT<yi2 — hTSi2) 
q{q + E„,-{M, + Mf)/2) 



0'12 = (7l • (72 , 5*12 — 3(Ti ■ q(J2 ■ q — CF 12 



Here hp = gy and 



9\ 



1 



3q 



1 

3 \q 



hp = gA 



1 



3q 



(53) 

For simplicity the dependence of the form factors gv and gA is not indicated and 
the terms containing l/rrip are omitted. The parts containing q^ + come from 
the induced pseudoscalar form factor gp for which the partially conserved axial- vector 
hypothesis (PCAC) is used. 

In order to calculate the nuclear matrix element in the coordinate space, one has to 
evaluate first the "neutrino potentials" that, at least in principle, depend explicitly on 
the energy E^^ of the virtual intermediate state 



(54) 



HK{r,2) = ^^Rj^ M^--), + ^„^_(M. + Af-/)/2 



where K = F,GT,T and fp^cT = io('i''^i2) and fp = —j2{q'fi2), and ri2 is the internu- 
cleon distance. 



With these "neutrino potentials" and the spin dependence given in eq. ( 52 1 where, 



naturally, in the tensor operator the unit vector q is replaced by fi2, the nuclear matrix 
element is now written (see 17^ ) 



(55) 



where |/) and \i) are the ground state wave functions of the final and, respectively, initial 
nuclei. The somewhat awkward definition of Af "^^ is used so that, if needed, one can use 
an "effective" value of the axial coupling constant gA but still use the tabulated values 
of the phase space integral G°'' that were evaluated with gA — I, 25. 
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► Define a valence space 

► Derive an effective interaction 

► Build and diagonalize the 
Hamiltonian matrix. 



Fig. 13. - Schematic illustration of the basic procedures in the nuclear shell model. 



6. Nuclear structure issues 



We will discuss now the procedures to evaluate the ground state wave functions of 
the initial and final nuclei |/) and \i) and the nuclear matrix element, eq. (551. There 
are two complementary methods to accomplish this task, the nuclear shell model (NSM) 
and the quasiparticle random phase approximation (QRPA). Since my own work deals 
with the QRPA method and the NSM will be covered by Prof. Poves in his seminar, the 
NSM will be described only superficially and reference to its results will be made mainly 
in comparison with the QRPA. 



6'1. Nuclear shell model. - The basic idea is schematically indicated in Fig. 13 The 
procedure should describe, at the same time, the energies and transition probabilities 
involving the low-lying nuclear states as well as the (3 and /3/3 decay nuclear matrix 
elements. The hamiltonian matrix is diagonalized in one of two possible bases, either 



(56) m - scheme = n„ij™^a+|0) = af^ ■ ■ ■ a+^|0) , D 




which is simple, based on the Slater determinants, and the corresponding hamiltonian 
matrix is sparse, but it has a huge dimension (d^, di, are the dimensions of the proton and 
neutron included subshells and p, n are the numbers of valence protons and neutrons) . 
The other possible basis uses states coupled to a given total angular momentum J and 
isospin T. The dimensions are smaller in that case but the evaluation of the hamilto- 
nian matrix element is more complicated and the corresponding hamiltonian matrix has 
fewer zero entries. Due to the high dimensionality of the problem one cannot include in 
the valence space too many single-particle orbits. Even the most advanced evaluations 
include just one oscillator shell, in most (3(3 decay candidate nuclei the valence space 
usually omits important spin-orbit partners. (For ^®Ge for example, the valence space 
consists of Pi/2,P3/2j /5/2 and gg/2 orbits for protons and neutrons, while omitting the 
essentially occupied /7/2 and empty 37/2 spin-orbit partners.) 
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Nucleus ^^"Sn ^^"5n "^Sb "-^ 7e "^Te 

Transition 0+ ^1+ 0+ ^ 1+ 4+ 3,4,5+ 0+ ^1+ |^ ^ |, |, 

Ti/jexp. 59.07m 3.72m 2.79m 3. 2d 12.5m 

Ti/2calc. (0.74) 32.21m 2.47m 1.56m 1.73d 6.42m 

Renorm. 0.54 0.6 0.55 0.54 0.53 



1347-^ 


i35Xe 


Cs 


0+ ^ 1+ 


3+1 3 5+ 
2 2 ' 2 ' 2 


5+ 4,5,6+ 


41.8m 


9.14h 


13.16d 


29.19m 


7.07h 


8. Id 


0.62 


0.63 


0.57 



Fig. 14. - Experimental and calculated /3-decay half-lives for several nuclei with A = 128-136 
(adopted from ref. [56]). 



In order to evaluate M^"^ the closure approximation is used and thus one evaluates 

(57) (/||O^IN>> with 0^'-E,,,,M/4f [(a+a+)^(a,aOlV 

where the creation operators create two protons and the annihilation operators annihilate 
two neutrons. In this way the problem is reduced to a standard nuclear structure problem 
where the many-body problem is reduced to the of evaluation of two-body transition 
densities. The matrix elements W^j'^ involve only the mean field (usually harmonic 
oscillator) one-body wave functions, and the transition operator defined above. 

The hamiltonian H^^f f is an effective operator. It is not entirely based on first principle 
reduction of the free nucleon-nucleon potential. Instead it relies on empirical data, in 
particular on the energies of states in semi-magic nuclei. At the same time effective 
operators should be used, in principle. Again, there is no well defined and well tested 
procedure to obtain their form. Instead, one uses empirical data and effective couplings 
(effective charges for the electromagnetic transitions, effective gA values for the weak 
transitions). 

How this procedure works is illustrated in Fig. [14] where the experimental f3 decay 
half- lives for a number of nuclei are compared with calculation and a quenching factor. 



i.e. the reduction factor of of ~0.57 is obtained. And in Fig. 15 we show the same 
comparison for the evaluation of the 2i//3P half-lives. (Note that this table is a bit obsolete 

1/2 



as far as the experimental data are concerned. In the meantime, the for ^-^"^Te was 



directly determined as 7.6x10^° years, and the limit for Tif2 in ^■^^Xe is longer than 
indicated, > 1-0 ^ 10^^ years, see the Table in the preceding section.) 

Various aspects of the application of the NSM to the /3/3 decay were reported in a 
number of publications [STJ [551 IMl ISHl EH [S3] • We will return to these results later 
when we compare the NSM and QRPA methods. 
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Parent nuclei 




76 Ge 


82 Se 


130 7-e 


I36xe 


Tlj^(g.s.) th. 
Tfj^ig.s.) exp 


3.7E19 


1.15E21 


3.4E19 


4E20 


6E20 


4.2E19 


1.4E21 


8.3E19 


2.7E21 


> 8.1E20 



Fig. 15. - Experimental and calculated 2v/3(3-dec&y half-lives for several nuclei (adopted from 
ref. [SS]). 



Here we wish to stress one general result of NSM evaluation of M^'^ . In the NSM one 
can classify states by their "seniority", i.e. by the number of valence nucleons that do 
not form Cooper- like pair and are therefore not coupled to /'^ = 0+. The dimensionality 
of the problem increases fast with seniority and thus it is of interest to see whether a 
truncation in seniority is possible and in general how the magnitude of M^'^ behaves as 



a function of the included seniority. This is illustrated in Fig. 16 from Ref. for the 
case of the ^^Ge decay. Note that the result does not saturate until — 12 which is 
nearly the maximum seniority possible. That behavior is observed in other cases as well 
(the saturation in Sm is faster in cases when one of the involved nuclei is semimagic) , 
and is relevant when the comparison to QRPA is made. 

NSM can also successfully describe nuclear deformation, unlike the usual application 
of QRPA. Many of the considered P(3 decay candidate nuclei are spherical or nearly so. 
However, one seemingly attractive candidate nucleus, ^^^Nd, is strongly deformed. More- 
over, the final nucleus, i^"Sm, is considerably less deformed than ^^^Nd. Unfortunately, 
the NSM is unable to describe this system, the dimension is just too large. However, 
to see qualitatively what the effect of deformation is, or the difference in deformation. 




Maximum seniority s^^^ 



Fig. 16. - The full matrix element M for the Ge — > Se transition evaluated with the 
indicated maximum seniority Sm included (adopted from ref. [63]). 
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• 1860=<Q>(Kr) 

• 2480 

• 2945 

• 3160 
U» 3300 

3620 
1010 
550 

^ I , L 

0.5 



A 



1 



Fig. 17. - The matrix element Af for the Se Kr transition for a large number of added 
quadrupole-quadrupole interaction strengths that induces increased deformation in *^Se. The 
magnitude of M"'' is ploted against the difference in the squared quadrupole moments of the 
invloved nuclei, (adopted from ref. [63] ) . 



might be, in Ref. [gSj the transition ^^Se ^^^Kr was considered and the initial nucleus 
®^Se was artificially deformed by including in the hamiltonian an additional (unrealis- 
tic) quadrupole-quadrupole interaction of varying strength. The result of that exercise 
are shown in Fig. [17] One can see that as the difference in deformation increases that 
magnitude of il/"'' decreases considerably. 

What are, then, the outstanding issues for NSM vis-a-vis the evaluation of the Af^" 
matrix elements? The most important one, in my opinion, is the limited size of the 
valence space. It is not clear how large (or small) the effect of the additional orbits 
might be. They cannot be, at present, included directly. Perhaps a perturbative method 
of including them can be developed, either in the shell model codes or in the definition 
of the effective operator. The other issue, that perhaps could be overcome, is the fact 
that the H^ff has not been determined for nuclei with A = 90 — 110 and thus the M^'^ 
are not available. Among them are important /3/9 candidate nuclei ^°°Mo and ^^Zr (a 
definitive calculations for ^^^Cd were not reported as yet either). 

6'2. QRPA basics. - The QRPA method was first applied to the charge changing 
modes by Halbleib and Sorensen [64, long time ago and generalized to include the particle- 
particle interaction by Cha [SS] . The use of quasiparticles in QRPA makes it possible to 
include the pairing correlations in the nuclear ground states in a simple fashion. With 
pairing included, the Fermi levels for protons and neutrons not only become diffuse, but 
the number of nucleons in each subshell will not have a sharp value, instead only a mean 
occupancy of each subshell will have a well determined value. 

To include the pairing we perform first the Bogoliubov transformation relating the 
particle creation and annihilation operators aj„,ajm with the quasiparticle creation and 
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-\~- UjCjui 



annihilation operators Cj^^^Cjm 
(58) 

where ajm — (— 1)-'^™%-™ and u'j 

The ampHtudes Uj,Vj are determined in the standard way by solving the BCS gap 
equations, separately for protons and neutrons, 

(59) = (2j„ + i)-i/2s,(2j, + i)i/v«c(j1;0+||F||j1;0+) , N^j:,{2j, + l)vl . 

Here N is the number of neutrons or protons and the gaps A are empirical quantities 
deduced from the usual mass differences of the corresponding even-even and odd-A nu- 
clei. We renormalize the strength of the pairing interaction (the coupling constant in 
Oa! ^'''II^IIJc; 0^)) shghtly such that the empirical gaps A are correctly reproduced. 

Our goal is to evaluate the transition amplitudes associated with charge changing 
one-body operator T"'*^ connecting the 0"*" BCS vacuum \0) of the quasiparticles c and 
in the even-even nucleus with any of the J'^ excited states in the neighboring odd-odd 
nuclei. In the spirit of RPA we describe such states as harmonic oscillations above the 
BCS vacuum. Thus 



\^QF!.Pa) 



(60) \J^M; m) = Sp„ [x'^^^j^A^pn; r M) + Y^^^j^Aipn; r M) 
where 

(61) A\pn- rU) = S,„^,„„ {jpfnp,j,,m^\JM)cl^^/^^^^ 

A{pn; ,rM) = {-ly-^' A{pn; .r - M) , 

and Oqupj^) is the "phonon vacuum" a correlated state that has the zero-point motion 
corresponding to the given built into it. It contains, in addition to the BCS vacuum, 
components with 4, 8, etc. quasiparticles. 

The so-called forward- and backward-going amplitudes X and Y as well as the corre- 
sponding energy eigenvalues are determined by solving the QRPA eigenvalue equa- 
tions for each J'^ 



(62) 



It is easy to see that eq. (62 1 is, in fact, an eigenvalue equation for of the type 
{A^ — B^)X — uP'X. Hence the physical solutions are such that uP' is positive, and we can 
choose uj to be positive as well. On the other hand, there could be unphysical solutions 
with < and hence imaginary energies. By varying the coupling constants in the 
matrices A and B we might trace the development of the solutions from the physical ones 
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with w > to the situation where w = 0. That point signals the onset of region where 
the original RPA (or QRPA) is no longer applicable, because the ground state must be 
rearranged. Examples are the transition from a spherical to deform shape, or transition 
from pairing of neutrons with neutrons (and protons with protons) to a pairing involving 
neutron-proton pairs. We will see that real nuclei are rather close to that latter situation, 
and we need to worry about the applicability of the method in such situations. 

To obtain the matrices A and B one needs first to rewrite the hamiltonian in the 
quasiparticle representation. Then 

(63) ^p^„,p,„, = (0|(ctct)(^*^)^iJ(4,4,)('^"'')|0> 

+ {UpVnUp'Vn' + VpUnVp'Un')gph{pn~'^ , J\V\p'n'~'^ , J) 

)9pp{pn,J\V\p'n',J) , 

and 

(64) B/„,p,„, ^ (0|i/(ctct)('^-*^)(-l)^^(ct,ci,)('^*'^|0) 

+ {~iy{UpVnVp'Un' + VpUnUp' Vn') dphipn^^ , J\V\p'n'^^, J) 
-{-iy{UpUnVp'Vn' + VpV.nUp>Un')gpp{pn, J\V\p'n', J) . 



Here Ep,En are the quasiparticle energies. The particle-hole and particle-particle 
interaction matrix elements are related to each other by the Pandya trasformation 

(65) {pn-\ J\V\yn'-\j) = -Ej,(2J' + 1) { j, } (P"', J\VWn, J) 



Above, in eqs. (64 65 1 we have introduced adjustable renormalization constants gph 
and gpp that multiply the whole block of interaction matrix elements for particle-hole 
and particle-particle configurations. Typically, one uses a realistic interaction for V (G- 
matrix) and thus the nominal values are gp^ = gpp = 1. 

The vectors X and Y obey the normalization and orthogonality conditions (for each 

•n 



(66) 



^pn^pn^pn 



pn pn 



^m^pn^p'n' '^pn^p'n' "pn,p'n' 



^pn^pn^ pn 

y^ Yni 
^m^pn^ p'n' 



T^m -v-m _ n 
pn pn ^ 

y-m y-m _ q 
pn p'n' ^ 



For a one-body charge-changing operator T'' the transition amplitude connecting 
the ground state of an even-even {N, Z) nucleus to the mth excited state in the odd-odd 
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{N -1,Z +1) nucleus is 
(67) (m; JM\T^'^^ 

where 



(68) 



y(2JTT) 



pn,j7r 



pn 



{-ly 



IT' 



The transition in the opposite direction, to the {N + 1, Z — 1) is governed by analogous 
formula but with X and Y interchanged. 

Given these formulae, we are able to calculate, within QRPA, the f3P decay nuclear 
matrix elements, for both 2j/ and modes. However, one needs to make another approx- 
imation, since the initial state jOg^p^; i) and the final state jOg^p^; /) are not identical. 
This is a "two vacua" problem. The standard way of accounting for the difference in the 
initial and final states is to add the overlap factor 



(69) {Jk\Jm) — ^pnXpnX^ 



pn pn ' 



where X, Y are the solutions of the QRPA equations of motion for the final nucleus. 
We can now write down the formula for the 2z^/3/3-decay matrix element as 



(70) 



m2^^^. (/II^^+II4>(i^IU.>(i+II'?^+IK> 



-+||i) = 'Spn{p\\a\\n){upVnX^„ + WpM„Fp") 
> = Sp„(p||a||n)(«pu„l^„ + up5„fl) • 



For the Ov decay the corresponding formula is 



(71) 



K 



J'^ ,ki ,kf ,^7" pnp' 7i' 



L Jn' Jp' J 

(p(l),p^);^ II /(ri2)r+T2+0,^/(ri2) II n(l), n'(2); J) x 

(o;il[c+i^.']j||J"M('/'^%l'/"^«>(^"fc/*ll[<£n]j||o+). 

The operators O/^, K — Fermi (F), Gamow- Teller (GT), and Tensor (T) contain neutrino 



The 



potentials, see eq. (54 1, and spin and isospin operators, and RPA energies 
reduced matrix elements of the one-body operators [CpC„]j in eq. (71 1 depend on the 
BCS coefficients u^, Vj and on the QRPA vectors X, Y\ they are just the reduced matrix 
elements of the one-body operators like in eqs. (67 68). The function j(rv2) in above 
represents the effect of short range correlations. These will be discussed in detail in 
the next Section. Note, that the radial matrix elements in eq. (71 1 are evaluated with 
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unsymmetrized two-particle wave functions. This is a generic requirement of RPA-like 
procedures as explained in Ref466 

6'3. Generalization - RQRPA. - The crucial simplifying point of QRPA (or RPA-like 
procedures in general) is the quasiboson approximation, the assumption that the com- 
mutation relations of a pair of fermion operators can be replaced by the boson commuta- 
tion relation. That is a good approximation for harmonic, small amplitude excitations. 
However, when the strength of the attractive particle-particle interaction increases, the 
number of quasiparticles in the correlated ground state I^qupa) increases and thus the 
method violates the Pauli principle. This, in turn, leads to an overestimate of ground 
state correlations, and too early onset of the QRPA collapse. 

To cure this problem, to some extent, a simple procedure, so-called renormalized 
QRPA (RQRPA) has been proposed and is widely used jH?]. In RQRPA the exact 
expectation value of a commutator of two-bifermion operators is replaced by 

{^qrpa\ [Mpn, JM),A+{p'n', JM)] |0+^p^) = V^™' x 



pn.J^ 



(72) 



with jp = y/2 jp + 1. 

To take this generalization into account, i.e. the nonvanishing values of Ppn.j" — 1 
one simply needs to use the amplitudes 

which are orthonormalized in the usual way instead of the standard X and Y everywhere 



also in the QRPA equation of motion ( 62 1 . In addition, the matrix elements of the one- 



1/2 

68|), must be multiplied by 1?^^ evaluated for the initial and 



body operators, eqs. (67 
final nuclei. 

1 /2 

In order to calculate the factors I?p„ one has to use an iterative procedure and 
evaluate in each iteration 

(74) = 1 - 2^^"'^P"' (^•/^fe(2-^ + ^Wpn'?) 

~ 2j^^P'^P''' (Sxfc(2J+ l)|Fptp) . 

Note the summation over the multipolarity J. Even if for the 21/ decay only J'^ = 1 + 
are seemingly needed, in RQRPA the equations need to be solved for all multipoles. 
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7. Numerical calculations in QRPA and RQRPA 

In the previous section all relevant expressions were given. But that is not all one 
needs in order to evaluate the nuclear matrix elements numerically. One has to decide, 
first of all, what are the relevant input parameters, and how to choose them. 

The first one to choose, as in all nuclear structure calculations, is the mean field 
potential, and which orbits are going to be included in the corresponding expressions. A 
typical choice is the calculation based on the Coulomb corrected Woods-Saxon potential. 
However, sometimes it is advisable to modify the single particle energy levels in order 
to better describe certain experimental data (energies in the odd-A nuclei or occupation 
numbers) . 

Next one needs to choose the nucleon-nucleon potential. For that one typically uses a 
G-matrix based on a realistic force. For example, in |66j the G-matrix used was derived 
from the Bonn-CD nucleon-nucleon force. It turns out that the results only weakly 
depend on which parametrization of the nucleon-nucleon interaction is used. 

Next, the effective coupling constants gph and gpp in the matrices A and B in equation 



of motion ( 62 1 must be determined. The parameter gph is not controversial. It is usually 
adjusted by requiring that the energy of some chosen collective states, often Gamow- 
Teller (GT) giant resonances, is correctly reproduced. It turns out that the calculated 
energy of the giant GT state is almost independent of the size of the single-particle basis 
and is well reproduced with gpii « 1. Hence, the nominal and unrenormalized value 
gph = 1 is used in most calculations; that reduces the number of adjustable parameters 
as well. 

The choice of the particle-particle parameter gpp is, however, not only important, but 
also to some extent controversial. One of the issues involved is illustrated in Fig. [18] 
There, one can see that the curves of M^'' versus gpp are rather different for different 
number of included single-paricle states. Thus, the calculated magnitude of AP^ will 
change dramatically if gpp is fixed and different number of levels is included. 



The other feature, illustrated in Fig. 18 is the crossing of M of zero value for 



certain gpp that is relatively close to unity. This was first recognized long time ago in 
Ref. [69^. Obviously, if AP" = then the 2i/ (3 f3-decay lifetime is infinite, this is a 
absolute suppression of that mode. Past that zero crossing the curves become very steep 
and the collapse of QRPA is reached when the slope becomes vertical. Obviously, the 
zero crossing is moved to larger values of gpp in RQRPA. 

One can use all of this to an advantage by abandoning the goal of predicting the M^'' 
values, but instead of using the experimental IVP'^ and determine the gpp (for a given set 
of s.p. states) in such a way, that the correct M^"^ is obtained. That is illustrated in Fig. 



19 The resulting gpp then depends on the number of included single-particle states, and 
is typically in the range 0.8 < gpp < 1.2 when a realistic G-matrix based hamiltonian is 
used. 

Adjusting the value of gpp such that the M'^" is correctly reproduced has been criti- 
cized, e.g. in Ref. |70j . There an alternative method, based on the experimentally known 
/3-decay ft values connecting the ground state of the intermediate nucleus (if that hap- 
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pn-QRPA pn-RQRPA 




Fig. 18. - The matrix elements M^" for ^'^Ge evaluated in QRPA (panel a) and RQRPA (panel 
b). The calculation was performed with the indicated number of included single-particle states. 
(Adopted from [SH]). The experimental value of M^"" is 0.13 MeV~^ 



pens to be 1+, which is so only in ^'^'^Tc, ^^^In and ^^^I among the /3/3-decay candidates). 
In QRPA the transition amplitudes for the EC process (decreasing nuclear charge) and 
(increasing nuclear charge) move in opposite way when gpp is increased; the first 
one goes up while the second one goes down. It is difficult, and essentially impossible 
to describe all three experimental quantities with the same value of gpp (in the case of 
lo^Mo this was noted already in [71 ). 

While the differences between these two approaches are not very large (see [T^ ) , there 




Fig. 19. - The calculated 2^ matrix elements for ^''Ge and *^Se as a function of gpp (solid lines). 
The experimental values are indicated by the horizontal dashed lines. The dotted line indicates 
the contribution of the first 1''' state. (Adopted from[70|.) 
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Fig. 20. - Running sum of the 2!//3/3-decay and Oi^/3/3-decay (only 1"*" component) matrix elements 
for ^^Ge and ^"''Mo (normalized to unity) as a function of the excitation energy Eex = En — 
{E, + Ef)/2. 



are other arguments why choosing the agreement with M^'^ is preferable. First, it is not 
really true that the first 1+ state is the only one responsible for the 2z//3/3 decay. This is 
illustrated for the cases of ^^Ge and ^'^"Mo in Fig. 20 Even though for ^"^^Mo the first 
state contributes substantially, higher lying states give non-negligible contribution. And 
in ^^Ge many 1+ states give comparable contribution. Thus, to give preference to the 
lowest state is not well justified, the sum is actually what matters. At the same time, the 
dilemma that the (3~ and (3~^ /EC matrix elements move with gpp in opposite directions 
makes it difficult to choose one of them. It seems better to use the sum of the products 
of the amplitudes, i.e. the 2i^(3(3 decay. 

At the same time, the contribution of the 1"*" multipole to the Qvj3f3 matrix element 
and the corresponding 2i'f3P matrix element are correlated, even though they are not 



identical, as also shown in Fig. 20 Making sure that the 2i//3/3 matrix element agrees 
with its experimental value constrains the 1^ part of the Qi'f3j3 matrix element as well. 

Since we are really interested in the M^^ matrix element, it is relevant to ask why 
do we fit the important parameter gpp to the 2v decay lifetime. To understand this, 
it is useful to point out that two separate multipole decompositions are built into Eq. 
(71 1. One is in terms the J'^ of the virtual states in the intermediate nucleus, the good 
quantum numbers of the QRPA and RQRPA. The other decomposition is based on the 
angular momenta and parities Sf'^ of the pairs of neutrons that are transformed into 
protons with the same J'^ . 

In Fig. 21 we show that it is essentially only the 1+ multipole that is responsible for 
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Fig. 21. - The contributions of difTerent intermediate-state angular momenta to M in 
^'^'^Mo (positive parities in the upper panel and negative parities in the lower one). We show 
the results for several values of Qpp. The contribution of the l"*" multipole changes rapidly with 
Qpp, while those of the other multipoles change slowly. 



the variation of Al'^'^ with gpp. (Note that in the three variants shown the parameter 
gpp changes only by 5%.) Fixing its contribution to a related observable {2v decay) 
involving the same initial and final nuclear states appears to be an optimal procedure for 
determining Qpp. Moreover, as was shown in |731 172] . this choice, in addition, essentially 
removes the dependence of M^^ on the number of the single-particle states (or oscillator 
shells) in the calculations. 

7'1. Competition between 'pairing' and 'broken pairs' . - The decomposition based on 
the angular momenta and parities J^^ of the pairs of neutrons that are transformed into 



protons with the same J'" is particularly revealing. In Fig. 22 we illustrate it both in the 
NSM and QRPA, with the same single-single particle spaces in each. These two rather 
different approaches agree in a semiquantitative way, but the NSM entries ior J > Q are 
systematically smaller in absolute value. There are two opposing tendencies in Fig. [22] 
The large positive contribution (essentialy the same in QRPA and NSM) is associated 
with the so-called pairing interaction of neutrons with neutrons and protons with protons. 
As the result of that interaction the nuclear ground state is mainly composed of Cooper- 
like pairs of neutrons and protons coupled to J — Q. The transformation of one neutron 




Fig. 22. - Contributions of different angular momenta J associated with the two decaying 
neutrons to the Gamow- Teller part of Af"" in *^Se (upper panel) and '^^"Te (lower panel). 
The results of NSM (dark histogram) and QRPA treatments (lighter histogram) are com- 
pared. Both calculations use the same single-particle spaces: (/5/2,P3/2, P1/2, 59/2) for *^Se and 
(57/2, ^5/2) ^3/2, •Si/2) ''■11/2) for ^^°Te. In the QRPA calculation the particle-particle interaction 
was adjusted to reproduce the experimental 2i//3/3-decay rate. 



Cooper pair into one Cooper proton pair is responsible for the J = ^ piece in Fig. |22[ 

However, the nuclear hamiltonian contains, in addition, important neutron-proton 
interaction. That interaction, primarily, causes presence in the nuclear ground state of 
"broken pairs", i.e. pairs of neutrons or protons coupled to J Their effect, as seen 
in Fig. 22 is to reduce drastically the magnitude of A'P^ . In treating these terms, the 
agreement between QRPA and NSM is only semi-quantitative. Since the pieces related 
to the "pairing" and "broken pairs" contribution ale almost of the same magnitude but of 
opposite signs, an error in one of these two competing tendencies is enhanced in the final 
M^'^ . The competition, illustrated in Fig. 22 is the main reason behind the spread of 
the published M"" calculations. Many authors use different, and sometimes inconsistent, 
treatment of the neutron-proton interaction. 

There are many evaluations of the matrix elements M^" in the literature (for the latest 
review see |21j). However, the resulting matrix elements often do not agree with each 
other as mentioned above and it is difficult, based on the published material, to decide 
who is right and who is wrong, and what is the theoretical uncertainty in Af*'"'. That 
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BKS-Ol = A. Bobyk, W. Kaminski, F. Simkovic, PRC63 (2001) <= 2yyS/3 20 times too slow 
Shell Model = E. Caurier, G. Martinez-Pinedo, F. Nowacki, A. Poves A. Zuker, RMP77 (2005) 
RFSV-07 = V.R., A. Faessler, R Simkovic, R Vogel, NPA793 (2007) (2003) <= 2y/3/3 fitted 
TFSG-86 = T. Tomoda, A. Faessler, K. W. Schmid, F. Grummer, NPA452 (1986) 

<= 2vJ3j3 8 times too fast 



Fig. 23. - Histogram of older published calculated values of {M^'^)^ for ''^Ge. The failure of 
some of the calculations to reproduce the known 2i//3/3-decay lifetime is indicated. 



was stressed in a powerful way in the paper by Bahcall et al. few years ago f74| where a 
histogram of 20 calculated values of (M"")^ for ^^Ge was plotted, with the implication 
that the width of that histogram is a measure of uncertainty. That is clearly not a valid 



conclusion as one could see in Fig. 23 where the failure of the outliers to reproduce the 
known 2i//3/3-decay lifetime is indicated. 

To see some additional reasons why different authors obtain in their calculations 
different nuclear matrix elements we need to analyze the dependence of the M^"^ on the 
distance r between the pair of initial neutrons (and, naturally, the pair of final protons) 
that are transformed in the decay process. That analysis reveals, at the same time, the 
various physics ingredients that must be included in the calculations so that realistic 
values of the M'^'^ can be obtained. 

7'2. Dependence on the radial distance. - The simplest and most important neutrino 



potential has the form (already defined in eq. (41 )) 



R 

(75) H(r) $(wr) , 

r 

where R is the nuclear radius introduced here as usual to make the potential, and the 
resulting M^"^, dimensionless (the 1/R^ in the phase space factor compensates for this). 
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Fig. 24. - The dependence on r of M°" for ^''Ge, ^'^^Mo and "°Te. The upper panel shows 
the full matrix element, and the lower panel shows separately 'pairing' [J' = of the two 
participating neutrons) and 'broken pair' 7^ 0) contributions. 



r is the distance between the transformed neutrons (or protons) and ^(wr) is a rather 
slowly varying function of its argument. 

From the form of the potential H(r) one would, naively, expect that the characteristic 
value of r is the typical distance between the nucleons in a nucleus, namely that r ^ R. 
However, that is not true as was demonstrated first in Ref. [66] and illustrated in Fig. 
[24] One can sec there that the competition between the "pairing" and "broken pairs" 
pieces essentially removes all effects of r > 2 — 3 fm. Only the relatively short distances 
contribute significantly; essentially only the nearest neighbor neutrons undergo the ov(3(i 
transition. The same result was obtained in the NSM [63]. (We have also shown in [66] 
that an analogous result is obtained in an exactly solvable, semirealistic model. There 
we also showed that this behaviour is restricted to an interval of the parameter that 
contains the realistic value near unity.) 

Once the r dependence displayed in Fig. |24]is accepted, several new physics effects 
clearly need to be considered. These are not nuclear structure issue per se, since they 
are related more to the structure of the nucleon. 

One of them is the short-range nuclcon-nucleon repulsion known from scattering ex- 
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without s.r.c. 
Spencer-Miler s.r.c. 
Co jastrow s.r.c. 
UCOM s.r.c. 




Fig. 25. - The r dependence of M^" in ^*^Ge. The four curves show the effects of different 
treatments of short-range correlations. The resulting M"" values are 5.32 when the effect is 
ignored, 5.01 when the UCOM transformation is applied and 4.14 when the treatment based 
on the Fermi hypernetted chain and 3.98 when the phenomenological Jastrow function is used. 
(See the text for details.) 



periments. Two nucleon strongly repel each other at distances r < 0.5 — 1.0 fni, i.e. the 
distances very relevant to evaluation of the M'^'^. The nuclear wave functions used in 
QRPA and NSM, products of the mean field single-nucleon wave function, do not take 
into account the influence of this repulsion that is irrelevant in most standard nuclear 
structure theory applications. The usual and simplest way to include the effect is to 
modify the radial dependence of the OvPf] operator so that the effect of short distances 
(small values of r) is reduced. This is achieved by introducing a phenomenological func- 
tion /(ri2 in the eq. (71 1. Examples of such Jastrow-like function were first derived in 
|75j and in a more modern form in |76j . That phenomenological procedure reduces the 
magnitude of M"^ by 20-25% as illustrated in Fig. |25) Recently, another procedure, 
based on the Unitary Correlation Operator Method (UCOM) has been proposed [77] , 
That procedure, still applied not fully consistently, reduces the M^" much less, only by 
about 5% |7S]. It is prudent to include these two possibilities as extremes and the corre- 
sponding range as systematic error. Once a consistent procedure is developed, consisting 
of deriving an effective Oj^/3/3 decay operator that includes (probably perturbatively) the 
effect of the high momentum (or short range) that component of the systematic error 
could be substantially reduced. 

Another effect that needs to be taken into account is the nucleon composite nature. 
At weak interaction reactions with higher momentum transfer the nucleon is less likely 
to remain nucleon; new particles for example pions, are often produced. That reduction 
is included, usually, by introducing the dipole form of the nucleon form factor, already 
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introduced, 



(76) fv.A 



1 + qyMt^.A 



where the cut- off parameters Mv,a have values (deduced in the reactions of free neutrinos 
with free or quasifree nucleons) ~ 1 GeV. This corresponds to the nuclon size of ~ 0.5—1.0 
fm. Note that in our case we are deaHng with neutrinos far off mass shell, and bound 
nucleons, hence it is not obvious that the above form factors are applicable. It turns out, 
however, that once the short range correlations are properly included (by either of the 
procedures discussed above) the M^" becomes essentially independent of the adopted 
values when Mv,a > 1 GeV. In the past various authors neglected the effect of short 
range correlations, and in that case a proper inclusion of nucleon form factor (or their 
neglect) again causes variations in the calculated M^" values. 

Yet another correction that various authors neglected must be included in a correct 
treatment. Since r < 2-3 fm is the relevant distances, the corresponding momentum 
transfer 1/r is of the order of ^200 MeV, much larger than in the ordinary /? decay. 
Hence the induced nucleon currents, in particular the pseudoscalar (since the neutrinos 
are far off mass shell) give noticeable contributions [75] . 

Finally, the issue of the axial current "quenching" should be considered. As shown in 



a rather typical example in Fig. 14 calculated values of the GT /3-decay matrix elements 
usually overestimate the corresponding experimental values. The reasons for that are, at 
least qualitatively understood, but the explanation is beyond the scope of these lectures. 
It suffice to say, that one can account for that effect phenomenologically and conveniently, 
by reducing the value of the axial-vector coupling constant gA to I instead of its true 
value 1.25. The phenomenon of quenching has been observed only in the GT /? decays, 
so it is not clear whether the same reduction of gA should be used also in the 0i//3f] 
decay. Nevertheless, it is prudent to include that possibility as a source of uncertainty. 



In anticipation , we already used the modified definition M"^" , see eq. ( 55 1 . 

We have, therefore, identified the various physics effects that ought to be included 
in a realistic evaluation of Af^" values. The spread of the calculated values, noted by 
Bahcall et al. [74] can be often attributed to the fact that various authors either neglect 
some of them, or include them inconsistently. 

8. Calculated M'"" values 

Even though we were able to explain, or eliminate, a substantial part of the spread 
of the calculated values of the nuclear matrix elements, sizeable systematic uncertainty 
remains. That uncertainty, within QRPA and RQRPA, as discussed in Refs. |731 [72] . is 
primarily related to the difference between these two procedures, to the size of the single- 
particle space included, whether the so-called quenching of the axial current coupling 
constant gA is included or not, and to the systematic error in the treatment of short 



range correlations In Fig. 26 the full ranges of the resulting matrix elements M"'' 
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Table II. - The calculated ranges of the nuclear matrix element M evaluated within both 
the QRPA and RQRPA and with both standard (qa ~ 1.254^ and quenched (qa ~ i-0) axtal- 
vector couplings. In each case we adjusted Qpp so that the rate of the 2uf3fi— decay is reproduced. 
Column 2 contains the ranges of M with the phenomenological Jastrow-type treatment of short 



range correlations (see I and II), while column 3 shows the UCOM-based results (see Ref. {ITp . 
Columns 3 and 5 give the Qu^jS-decay half-life ranges corresponding to the matrix- element ranges 



in columns 2 an 


d 4, for < 171(30 > 
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Nuclear 


(R)QRPA (Jastrow s.r.c.) 


(R)QRPA (UCOM s.r.c.) 


transition 


A/'"" rj'/a {(miifi) = 50 meV) 




rr/a iimp^) = 50 meV) 
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(3.33,4.68) 
(2.82,4.17) 
(1.01,1.34) 
(2.22,3.53) 
(1.83,2.93) 
(2.46,3.77) 
(2.27,3.38) 
(1.17,2.22) 


(6.01,11.9) X 10^*^ 
(1.71,3.73) X 10^*^ 
(7.90, 13.9) X 10^^ 
(1.46,3.70) X 10^** 
(1.95,5.01) X 10^'' 
(3.33,7.81) X 10^'' 
(1.65,3.66) X lO^'' 
(3.59, 12.9) X 10^'" 


(3.92,5.73) 
(3.35,5.09) 
(1.31,1.79) 
(2.77,4.58) 
(2.18,3.54) 
(3.06,4.76) 
(2.84,4.26) 
(1.49,2.76) 


(4.01,8.57) X 10^*^ 
(1.14,2.64) X 10^*^ 
(4.43,8.27) X 10^^ 
(8.69,23.8) X 10^^^ 
(1.34,3.53) X 10^*^ 
(2.09,5.05) X lO^'^ 
(1.04,2.34) X 10^*^ 
(2.32, 7.96) X 10^'^ 



is indicated. The indicated error bars are highly correlated; e.g., if true values are near 
the lower end in one nucleus, they would be near the lower ends in all indicated nuclei. 

The figure also shows the most recent NSM results [HS]. Those results, obtained with 
Jastrow type short range correlation corrections, are noticeably lower than the QRPA 
values. That difference is particularly acute in the lighter nuclei ''^Ge and ^^Se. While 
the QRPA and NSM agree on many aspects of the problem, in particular on the role 
of the competition between "pairing" and "broken pairs" contributions and on the r 
dependence of the matrix elements, the disagreement in the actual values remains to be 
explained. 

When one compares the and Oi/ matrix elements (Figs. 10 and 26) the feature to 
notice is the fast variation in M^'' when going from one nucleus to another while M°'^ 
change only rather smoothly, in both QRPA and NSM. This is presumably related to the 
high momentum transfer (or short range) involved in Qv[3(3. That property of the M*^" 
matrix elements makes the comparison of results obtained in different nuclei easier and 
more reliable. 

Given the interest in the subject, we show the range of predicted half-lives corre- 
sponding to our full range of M'^" in Table [n] (for {mpp) = 50 mcV). As we argued 
above, this is a rather conservative range within the QRPA and its related frameworks. 
One should keep in mind, however, the discrepancy between the QRPA and NSM results 
as well as systematic effects that might elude either or both calculations. Thus, as we 
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Fig. 26. - The full ranges of M "'^ with the two alternative treatments of the short range 
correlations. For comparison the results of a recent Large Scale Shell Model evaluation of AI ""^ 
that used the Jastrow-type treatment of short range correlations are also shown (triangles). 



have seen, while a substantial progress has been achieved, we are still somewhat far from 
being able to evaluate the nuclear matrix elements confidently and accurately. 



Appendix 



Neutrino magnetic moment and the distinction between Dirac and Majorana neutrinos 
The topic of neutrino magnetic moment is seemingly unrelated to the Oi//9/3 decay. 
Yet, as will be shown below, experimental observation of allows one to make important 
conclusions about the Dirac versus Majorana nature of neutrinos. Hence, it is worthwhile 
to discuss it here. 

Neutrino mass and magnetic moments are intimately related. In the orthodox Stan- 
dard Model neutrinos have a vanishing mass and magnetic moments vanish as well. 
However, in the minimally extended SM containing gauge-singlct right-handed neutrinos 
the magnetic moment /x^ is nonvanishing, but unobservably small [80j . 

(A.l) Ai,= ^^^m, = 3x 



^87r2 '^"leV 



An experimental observation of a magnetic moment larger than that given in eq.(A.ll 
would be an uneqivocal indication of physics beyond the minimally extended Standard 
Model. 

Laboratory searches for neutrino magnetic moments are typically based on the ob- 
sevation of the v — e scattering. Nonvanishing will be recognizable only if the cor- 
responding electromagnetic scattering cross section is at least comparable to the well 
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understood weak interaction cross section. The magnitude of (diagonal in flavor or 
transitional) which can be probed in this way is then given by 

(A.2) ^ ^ ^ 10-10 

where T is the electron recoil kinetic energy. Considering realistic values of T, it would 
be difficult to reach sensitivities below ~ lO^^^/is using the v — e scattering. Present 
limits are about an order of magnitude larger than that. 

Limits on fj,i, can also be obtained from bounds on the unobserved energy loss in 
astrophysical objects. For sufficiently large Hi, the rate of plasmon decay into the vD 
pairs would conflict with such bounds. Since plasmons can also decay weakly into the 
pairs , the sensitivity of this probe is again limited by the size of the weak rate, leading 
to 

(A.3) — = —= — ncjp , 

where wp is the p lasmon frequency. Since (hujp)'^ ^ rn^T that limit is stronger than that 
given in eq.(A.2l. Current limits on /ij^ based on such considerations are ~ 10~^^^b- 

The interest in and its relation to neutrino mass dates from ~1990 when it was 
suggested that the chlorine datajlj on solar neutrinos show an anticorrelation between 
the neutrino flux and the solar activity characterized by the number of sunspots. A 
possible explanation was suggested in Ref. [81] where it was proposed that a magnetic 
moment /ij/ ^ 10~'-^°~^^V-B would cause a precession in solar magnetic field of the 
neutrinos emitted initially as left-handed Ve into unobservable right-handed ones. Even 
though later analyses showed that the effect does not exist, the possibility of a relatively 
large fXi, accompanied by a small mass rrii, was widely discussed and various models 
accomplishing that were suggested. 

If a magnetic moment is generated by physics beyond the Standard Model (SM) at 
an energy scale A, we can generically express its value as 

(A.4) ^"^1^' 

where e is the electric charge and G contains a combination of coupling constants and 
loop factors. Removing the photon from the diagram gives a contribution to the neutrino 
mass of order 

(A.5) ~ GA. 

We thus have the relationship 



A^ /i^ 

2me /is 10-18^^ I 



(A.6) -.-^^ - T7^[A(TeV)F eV, 



which implies that it is difficult to simultaneously reconcile a small neutrino mass and a 
large magnetic moment. These considerations are schematically illustrated in Fig. 27 
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Fig. 27. - a) Generic contribution to the neutrino magnetic moment induced by physics beyond 
the standard model, b) Corresponding contribution to the neutrino mass. The soHd and wavy 
lines correspond to neutrinos and photons respectively, while the shaded circle denotes physics 
beyond the SM. 



This nai've restriction given in Eq.(A.6l can be overcome via a careful choice for the 
new physics, e.g., by requiring certain additional symmetries [551 IMJ Ull US] • Note, how- 
ever, that these symmetries are typically broken by Standard Model interactions. For 
Dirac neutrinos such symmetry (under which the left-handed neutrino and antineutrino 
V and v"^ transform as a doublet) is violated by SM gauge interactions. For Majorana 
neutrinos analogous symmetries are not broken by SM gauge interactions, but are instead 
violated by SM Yukawa interactions, provided that the charged lepton masses are gen- 
erated via the standard mechanism through Yukawa couplings to the SM Higgs boson. 
This suggests that the relation between and is different for Dirac and Majorana 
neutrinos. This distinction can be, at least in principle, exploited experimentally, as 
shown below. 

Earlier, I have quoted the Ref.[52] (see Fig 4]) to stress that observation of the Qv(i(3 
decay would necessarily imply the existence oi a novanishing neutrino Majorana mass. 
Analogous considerations can be applied in this case. By calculating neutrino magnetic 
moment contributions to mi, generated by SM radiative corrections, one may obtain in 
this way general, "naturalness" upper limits on the size of neutrino magnetic moments 
by exploiting the experimental upper limits on the neutrino mass. 

In the case of Dirac neutrinos, a magnetic moment term will generically induce a 
radiative correction to the neutrino mass of order [85] 

(A.7) m,.-^ — ^ ^ — ^^^[A(TeV)f eV. 



Taking A ~ 1 TeV and m,^ < 0.3 eV, we obtain the limit /i^ < 10~^^/iB (and a more 
stringent one for larger A), which is several orders of magnitude more constraining than 
current experimental upper limits on /Xj,. 

The case of Majorana neutrinos is more subtle, due to the relative flavor symmetries 
of and /i^ respectively. For Majorana neutrinos the transition magnetic moments 
[Mi'Iq/s antisymmetric in the flavor indices {a,/3}, while the mass terms [toj/Jq/j are 
symmetric. These different flavor symmetries play an important role in the limits, and 
are the origin of the difference between the magnetic moment constraints for Dirac and 
Majorana neutrinos. 

It has been shown in Ref. [57] that the constraints on Majorana neutrinos are sig- 
nificantly weaker than those for Dirac neutrinos |86). as the different flavor symmetries 
of TOi/ and 111, lead to a mass term which is suppressed only by charged lepton masses. 
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Fig. 28. - One loop diagram contributions to the Majorana neutrino mass associated with the 
magnetic moment that sum to zero (see ^7])- The cross indicates the magnetic moment [ii\ap 
and the / is the lepton doublet. 



This conclusion was reached by considering one-loop mixing of the magnetic moment 
and mass operators generated by Standard Model interactions. The authors of Ref.^87] 
found that if a magnetic moment arises through a coupling of the neutrinos to the neutral 
component of the SU{2)i^ gauge boson, the constraints for /Xi-e and fj-rfj. are comparable 
to present experiment limits, while the constraint on /ig^t is significantly weaker. Thus, 
the analysis of Ref . |,87j lead to a model independent bound for the transition magnetic 
moment of Majorana neutrinos that is less stringent than present experimental limits. 



Those considerations are illustrated in Figs. 28 29 

Even more generally it was shown in Ref. [88 that two- loop matching of mass and 
magnetic moment operators implies stronger constraints than those obtained in[87 if the 
scale of the new physics A > 10 TeV. Moreover, these constraints apply to a magnetic 
moment generated by either the hypercharge or SU (2) ^ gauge boson. In arriving at 
these conclusions, the most general set of operators that contribute at lowest order to 
the mass and magnetic moments of Majorana neutrinos was constructed, and model 
independent constraints which link the two were obtained. Thus the results of Ref. |88j 
imply completely model independent naturalness bound that - for A > 100 TeV - is 
stronger than present experimental limits (even for the weakest constrained element 
/j,e^). On the other hand, for sufficiently low values of the scale A the known small values 
of the neutrino masses do not constrain the magnitude of the magnetic moment /i^ more 
than the present experimental limit. Thus, if these conditions are fulfilled, the discovery 
of /Liy might be forthcoming any day. 

The above result means that an experimental discovery of a magnetic moment near 
the present limits would signify that (i) neutrinos are Majorana fermions and (ii) new 





Fig. 29. - The one and two loop contributions to the Majorana neutrino mass associated with 
the magnetic moment. Here X is the charged lepton mass insertion (see [87]). 
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lepton number violating physics responsible for the generation of fii, arises at a scale A 
which is well below the see-saw scale. This would have, among other things, implications 
for the mechanism of the neutrinoless double beta decay and lepton flavor violation as 
discussed above and in Ref.^28j. 

* * * 

The original results reported here were obtained in collaboration with Nicole Bell, Vin- 
cenzo Cirigliano, Jonathan Engel, Amand Faessler, Michail Gorchtein, Andryi Kurylov, 
Gary Prezeau, Michael Ramsey-Musolf, Vadim Rodin, Fedor Simkovic, Peng Wang and 
Mark Wise. The fruitful collaboration with them is gratefully acknowledged. 
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